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NOTATION 


tip  loss  factor;  blade  elements  outboard  of  radius  BR 
are  assumed  to  have  profile  drag  but  no  lift. 

profile-drag  power  coefficient,  P  /[pwR* (flR) 3] 

o 

thrust  coefficient,  T/[pirRa  (QR)  a] 
parasite  drag  of  helicopter  fuselage 
endurance 

mass  moment  of  inertia  of  one  blade  about  a  flapping 
hinge  (1/3  mRa) 

lift 

Mach  number  at  radial  position  x,  azimuth  t|> 

power  available  at  main  rotor  shaft 

climb  power,  i.e.  excess  power  available  for  climb 

power  component  arising  from  compressibility  effects  on 
advancing  blade 

power  component  arising  from  rotor  lift  (Induced  Power) 

installed  power  i.e.  total  power  output  of  helicopter 
engines 

power  component  arising  from  rotor  drag  (Profile  Power) 

power  component  arising  from  parasite  drag 

power  required  at  the  main  rotor  shaft 

P.  +  P  +  P 
1  o  p 

blade  tip  radius 
specific  fuel  consumption 
main  rotor  thrust 

main  rotor  thrust  out  of  ground  effect 

resultant  velocity  perpendicular  to  blade-span  axis 
at  blade  element 

component  at  blade  element  of  resultant  velocity 
perpendicular  both  to  blade-span  axis  and  U 


NOTATION  (CONT.D) 


component  at  blade  element  of  resultant  velocity 
perpendicular  to  blade-span  axis  and  to  ANF 

airspeed  of  helicopter  along  flight  path 

airspeed  increment 

maximum  value  of  V 

rate  of  climb 

V/v. 

hov 

all  up  weight  of  helicopter 
weight  of  fuel 

slope  of  curve  of  section  lift-coefficient  against 
section  angle  of  attack 

constant  term  in  Fourier  series  that  expresses  @  hence, 
the  rotor  coning  angle 

coefficient  of  cosifi  in  Fourier  series  that  expresses  0 
hence,  first-order  longitudinal  flapping  angle 

speed  of  sound 

number  of  main  rotor  blades 

coefficient  of  sin4>  in  Fourier  series  that  expresses  6 
hence,  first  order  lateral  flapping  angle 

blade  section  cord 

equivalent  blade  chord,  (/R  cra  dr)/(/R  ra  dr) 

o  o 

blade  section  lift  coefficient 

average  section  lift  coefficient  in  reversed-velocity  region 

blade  section  profile-drag  coefficient 

average  section  profile-drag  coefficient  in  reversed 
velocity  region 

equivalent  flat-plate  area  of  helicopter 

semi-empirical  factor  for  induced  velocity  calculation  (Eq.9) 
mass  of  single  rotor  blade 


blade  element  radius 


NOTATION  (CONT.D) 


non-dimensional  resultant  velocity  at  blade  element,  U/QR 

non-dimensional  component  of  resultant  velocity  at  blade 

element,  U  /OR 
P 

non-dimensional  component  of  resultant  velocity  at  blade 
element,  U^/QR 

ratio  of  blade  element  radius  to  blade  tip  radius,  i.e.  r/R 
height  of  rotor  above  ground  level 

p/,pISA+15°C 

angular  velocity  of  rotor 

angle  of  attack  of  rotor,  i.e.  angle  between  axis  of  no 
feathering  and  plane  perpendicular  to  flight  path. 

(Positive  when  axis  is  pointing  rearward) . 

angle  of  attack  of  blade-element,  measured  from  line  of 
zero  lift 


ar  at  radial  position  x,  azimuth  ip 

blade  flapping  angle  at  particular  azimuth  position, 
measured  from  plane  perpendicular  to  ANF 

Lock's  number,  cpa  rVi 


climb  angle,  i.e.  angle  between  flight  path  and  horizontal 
(positive  when  climbing) 


coefficients  in  power  series  expressing  c.  as  a  function 

a 


of  a  (c. 
r  d 


50  +  6l°r  +  62°ra> 


o 

blade  section  pitch  angle,  i.e.  angle  between  zero  lift 
line  of  blade  section  and  plane  perpendicular  to  the  ANF 

blade  pitch  angle  at  hub 

blade  pitch  angle  at  tip 

blade  twist,  i.e.  9„,_  -  0 
TIP  o 

inflow  ratio,  (Vsina  -  v)/(QR) 

tip  speed  or  advance  ratio,  (Vcosa) / (DR) 

induced  flow  velocity  at  rotor  (always  positive) 


NOTATION  (CONT.D) 


momentum  theory  value  of  v  at  hover,  QR/C^/2 


factored  induced  velocity  i.e.  v  multiplied  by  a  semi- 
empirical  factor  to  allow  for  non-uniform  distribution 
of  v  over  the  rotor  disc 


density  of  air 

density  of  air  at  sea-level,  ISA+15°C  conditions 
ISA+15°C 

blade  solidity,  i.e.  the  ratio  between  blade  area  and 

disc  area,  be  /it R 
e 

inflow  angle  at  blade  element  in  plane  perpendicular  to 
blade-span  axis 

blade  azimuth  angle  measured  in  direction  of  rotation 
(zero  at  downwind  position) 


GLOSSARY 


Axis  of  No-Feathering;  defined  as  the  axis  about 
which  there  is  no  cyclic  pitch  change,  i.e.  an  observer 
rotating  on  this  axis  would  observe  blade  flapping 
but  no  pitch  change. 

Aircraft  Research  and  Development  Unit. 

International  Civil  Aviation  Organization. 

In  Ground  Effect. 

International  Standard  Atmosphere. 

Maneuver  Criteria  Evaluation  Program  (Reference  3  and  4 ) . 
No  Feathering  Plane;  plane  perpendicular  to  ANF. 

Out  of  Ground  Effect. 

Equivalent  sea-level  ambient  pressure  determined  from 
ambient  pressure  at  airfield  and  airfield  altitude. 

Royal  Australian  Navy. 

Systeme  Internationale  (International  System  of  units) . 
Specific  Fuel  Consumption  (applied  to  helicopter  engines 
in  this  case) . 

Tip  Path  Plane;  defined  as  the  plane  which  is  traced  out 
by  the  blade  tips  of  the  rotor,  i.e.  an  observer  on  this 
plane  would  observe  blades  changing  pitch  but  not  flapping. 


NOTE  ON  UNIT  SYSTEMS 


The  programs  described  are  capable  of  processing  data 
in  both  S.I.  and  Imperial  systems.  The  latter  is  the  major  system 
of  this  report  since  aeronautical  practice  still  makes  use  of  the 
Imperial  system.  The  programs  always  input  and  output  airspeed 
in  knots  and  power  in  horsepower  regardless  of  which  system  is 


1. 


INTRODUCTION 


In  early  1983,  the  Aircraft  Behaviour  Studies  -  Rotary  Wing 
Group  at  ARL,  was  asked  to  supply  performance  estimates  for  various 
helicopters  being  considered  for  purchase  by  the  RAN.  These  estimates 
included  power  required  in  steady  level  and  steady  climbing  flight, 
range,  endurance  and  ceiling  approximations.  Two  methods  of  performance 
estimation  were  readily  available: 

a)  Performance  charts  (Refs.l  and  2). 

b)  Program  'MCEP'  (Maneuver  Criteria  Evaluation 
Program,  Refs. 3  and  4). 

These  performance  charts ,  while  giving  an  indication  of  the 
onset  of  blade  stall,  have  some  limitations: 

(i)  Each  set  of  charts  apply  to  one  particular  aerofoil 
section;  Ref.l  uses  NACA  0012,  Ref. 2  uses  NACA  23012. 

(ii)  The  charts  are  restricted  to  certain  flight  conditions, 
which  may  necessitate  interpolation  between  a  series  of 
charts . 

(iii)  The  charts  are  sometimes  difficult  to  read  accurately, 
with  small  reading  errors  often  resulting  in  large 
discrepancies  in  the  power  components. 

(iv)  Many  of  the  charts  of  reference  1  apply  only  to  very 
high  performance  helicopters  (high  advance  ratios) 
which,  at  present,  do  not  exist. 

(v)  The  analysis  of  reference  2  uses  momentum  theory, 

which  assumes  uniform  inflow,  to  find  v.  In  practice, 
at  low  speed,  the  average  v  over  the  disc  is  higher 
due  to  the  effects  of  blade  tip-loss  and  non-uniform 
inflow.  Thus  the  charts  of  reference  2  underestimate 
the  induced  power. 

(vi)  For  some  performance  calculations  it  is  necessary  to 
produce  a  speed-power  polar.  Using  charts  to  create 
such  a  polar  would  be  a  slow  and  laborious  task. 


Program  'MCEP'  is  a  large  and  complex  program  which  can  be 
used  to  evaluate  a  wide  range  of  helicopter  maneouvres.  'MCEP'  can 
also  estimate  a  level  flight  speed-power  polar,  but  the  method  of 
outputting  data  is  unsatisfactory  for  present  applications.  For  this 
reason  program  'ENERGY'  has  been  produced  at  ARL  (see  Appendix  4) , 
which  uses  the  same  relatively  simple  expressions  for  power  estimation 
and  induced  velocity  as  'MCEP',  but  is  much  smaller  and  outputs  data 
in  a  more  suitable  form.  While  'ENERGY'  gives  accurate  results  in 
many  cases  (see  Ref. 4),  certain  deficiencies  are  evident. 

(i)  The  occurrence  of  blade  stall  is  not  recognized; 

(ii)  Blade  twist  is  not  taken  into  account. 

Power  variations  due  to  these  factors  would  therefore  not  be  predicted 
by  'ENERGY',  thus  causing  errors  in  performance  estimation. 

In  order  to  overcome  the  above-mentioned  deficiencies  of  these 
performance  estimation  methods,  program  'POLAR'  has  been  created,  which 
uses  a  simplified  blade-element  analysis,  based  on  the  work  of  Gessow  and 
Crim  (Ref.  5),  and  is  applicable  to  any  conventional  helicopter  (i.e. 
single  main  rotor)  with  flapping  blades. 

Due  to  assumptions  made  and  the  resulting  simplification  of 
the  analysis,  'POLAR'  is  applicable  for  uSQ.3.  It  should  be  noted 
that  this  limit  is  not  exceeded  by  most  conventional  helicopters 
operational  at  present. 

A  description  of  the  equations  used  and  assumptions  made 
is  given  in  Section  2,  while  Section  3  gives  a  description  of  'POLAR's 
program  structure.  Section  4  is  a  User's  guide  to  'POLAR'  which 
demonstrates  the  program's  use  in  estimating  helicopter  performance. 

A  comparison  between  'POLAR'  and  other' methods  of  performance  estimation 
is  provided  in  Section  5. 

2.  EQUATIONS 

This  section  presents  the  equations  used  by  program  'POLAR'. 

The  method  by  which  these  equations  are  solved  is  given  in  Section  3. 


[3] 


2.1  Forces  In  Steady  Flight 

Figure  1  shows  the  forces  acting  on  the  helicopter  in  steady 
powered  flight.  The  resultant  thrust  of  the  main  rotor  is  assumed  to 
be  perpendicular  to  the  Tip  Path  Plane.  Bramwell  (Ref. 6)  shows  that 
this  assumption  is  accurate. 

Resolving  forces  vertically  and  horizontally  gives*: 


W  +  Dp  sinyc  =  T  cos  [  (-a)  -  yc  "  a  ] 


Dp  cosy  =  T  sin  (  (-a)  -  y  -  a  I 

C  Cl 


From  equation  (1) , 


..(2) 


where 


W  +  Dp  sinYc 

cos  (a+Y  +a, ) 
c  1 


Dp  =  -  pV*f 


Equivalent  flat  plate  area,  f,  is  assumed  to  be  independent 
of  fuselage  angle  of  attack. 

Dividing  equation  (2)  by  equation  (1) ,  and  re-arranging 


From  momentum  theory: 


X  *  p  tan  a - - - ^ 

2 (pa  +  X*)4 

which ,  for  hover  (i.e.,p~0)  reduces  to 


and  in  forward  flight,  assuming  a  is  small  and  X<<p,  gives 
X  -  -  CT  /  (2p) 

After  the  initial  estimate  for  X  is  made  using  the  above 
equations,  v  (and  then  X,  see  §3.2)  is  calculated  according  to  the 
standard  resultant  velocity  equation: 

v  =  (  V  +  v  -  2  V  \J  sina  )  * 

where  v  -  v/v^ov,  V  *  V/\»hQv  and  vhQv  =  QR/CT/2  which  can  be  solved 
iteratively,  assuming  V  and  a  are  known. 

At  low  speeds  however,  the  value  of  v  obtained  above  is  not 
realised  in  practice,  because  the  induced  velocity  distribution  is 
non-uniform  over  the  rotor  disc.  Below  p  *  0.14,  a  semi -empirical 
factor  (Refs.  3  and  4)  is  thus  applied  to  the  value  of  v  to  give  a 
more  realistic  result. 

vf  ■  k.v 

where  k  *  1.0  for  p  >  0.14 

and  k«  1.0  -  2.14  (p  -  0.14)  for  p  &  0.14 

2.3  Blade  Flapping  Motion 

In  steady  flight,  the  blade-motion  must  be  periodic  and  is 
therefore  capable  of  being  expressed  in  a  Fourier  series.  Thus  for 
rigid  blades: 


[5] 


g  *  a^  -  costji  -  b^  siniji 

ignoring  terms  of  second  order  and  above. 

Thus  aQ  represents  the  blade  coning  angle, 

a^  represents  the  first  order  longitudinal  flapping  angle, 
b^  represents  the  first  order  lateral  flapping  angle. 

Reference  6  states  that  the  higher  harmonics  of  the  Fourier 

series  (a^,  a3'  etc.)  represent  blade  displacements  which  are  of 

the  same  order  as  elastic  deflections  of  the  rotor  blades.  Thus  it  is 
inconsistent  to  calculate  the  higher  harmonics  of  rigid  blade  motion 
without  including  structural  deflections  of  the  blade. 

2.4  Blade  Element  Angle  of  Attack 

Figure  2  represents  a  typical  blade  element.  Reference  7  shows 
that  for  small  g,  the  velocity  components  at  the  rotor  blade  element 
may  be  expressed  non-dimensionally  as: 

°T 

UT  *  OR  31  X  +  U  Sin'!' 

u  1 

up  *  OR  “  x  +  I  wai+  (_ya0  +  3d5l)  C°S* 

+  (-xa^)sini|/  +  j  ya1cos2iJ>  +  j  ub1sin2t|) 


Also  from  figure  2, 

<!>  a  arctan  Up 
UT 

Thus  the  local  angle  of  attack  of  a  blade  element  can  be 
found  from: 


<*r  "  9o  +  x0l  +  ♦ 


.. (10) 


K32Z 


itities  Derived  From  Blade  Elei 


The  extended  analysis  of  Gessow  and  Crim  (Ref.  5)  has 
been  used  as  the  basis  for  all  blade  element  equations  presented 
here.  This  extended  theory  differs  from  standard  rotor  theory  in 
that  no  limitation  is  placed  on  the  inflow  angle  and  the  reversed 
velocity  region  is  taken  into  account.*  However,  no  allowances 
are  made  for  blade  stall  outside  of  the  reversed  velocity  region, 
and  compressibility  effects  are  ignored. 

The  equations  presented  here  differ  from  those  used  in 
reference  5  in  that: 


(a)  Second  order  flapping  coefficients  are  neglected, 
for  reasons  stated  earlier. 

(b)  Terms  of  high  order  in  y  and  9^  cure  neglected. 

High  order  terms  in  y  are  ignored  since  above 
approximately  y  *  0.3,  blade  stall  and  compressibility 
effects  become  significant.  As  most  blades  in  present 
day  conventional  helicopters  have  twist  less  than 

12  deg,  high  order  terms  in  9^  are  ignored. 

These  modifications  vastly  simplify  the  equations  and 
method  of  solution.  By  ignoring  second  order  flapping,  for  example, 
the  five  simultaneous  flapping  equations  of  Ref.  5  reduce  to 
three  equations  which  can  be  solved  directly. 


ir  - 

I 

*  When  a  helicopter  is  in  forward  flight,  part  of  the  inner  section 
of  the  retreating  blade  is  rotating  in  the  aft  direction  at  a 
speed  which  is  lower  them  the  forward  speed.  Thus  air  flow  is  in 
a  reverse  direction  in  this  section,  which  is  known  as  the  reversed 
velocity  region. 


After  likewise  simplif ication ,  the  profile  power  coefficient 
is  given  by: 


0  <0  -  *1  (l  +  *  **(?♦!  "!)+  K3  (j*  f  -5) 

*‘4  0  +  B  »)♦(»  •  5  -DO  *5  *») 

+  I“2  1  K7  +  X!  (I  K19  *  T  K2o) 

*(K  +  T  bD  (4  K19  +  ?  K2o) 

where  K  =  5Q  +  ^  sin0Q  +  $2  sin20Q 


K2  *  K7  91 


K3  “  2  Ka  91 


K4  "  3  K9  91 


K5  "  K6  "  ° 


K?  =  cosQq  +  6 2  sin20Q 


Kg  *  (-  sin0Q  +  2  62  cos29q) 


k9  “  (“  |  \  cos90  '  2  62  sin290)  91 


K10  “  K11  *  K12  "  ° 


K13  “  52  COS  9C 


K,  ■  -  5„  sin  (29.)  0, 


9 


K19  “  60  +  51  Sin60  +  S2  COs2e0 


K20  -  (51  cos0o  -  <$2  sin290) 


K_,  «•  K__  ■  =  K„.  -  0 

21  22  23  24 


Certain  'K*  constants  which  are  equated  to  zero  result 
from  neglect  of  high  order  terms  in  0^  and  thus  have  been  ignored. 

and  K^4  are  used  in  Ref.  5  to  calculate  the  decelerating  torque. 
Since  accelerating  and  decelerating  torque  are  not  required  for 
performance  calculations  of  helicopters  in  power-on  flight,  they 
sure  not  considered  here. 


2 . 6  Power  Calculations 


The  power  absorbed  by  a  helicopter  is  made  up  of  several 
components .  The  standard  equations  used  to  calculate  these  sure  shown 


below. 


P.  ■  T  .  v_ 
l  f 


(Induced) 


P  ■  -  p  V3  f 
P  2 


(Parasite) 


P  *  C  (pwR* (HR) 3  ]  (Profile) 

°  *o 


where  Cp  is  given  by  equation  (15) . 
o 

For  convenience,  parameter  P  is  defined  as 

sum 


P  »  P.  +  P  +  P 
sum  i  p  o 


Climb  power  Pc  is  calculated  according  to  whether 


(a)  Climb  angle  y  is  known,  in  which  case  Ref.  2 

c 


gives  the  following  expression  - 


[10] 


p 

c 


sm> 


P 

-*lnYcife  » 


cosa 


l-cosaY 


V 


c  Ta(0R)a 


cosaa 

Ma 


*1 


_u _ 

cosa 


Conversely  the  climb  angle  may  be  determined  if  the  available  climb 
power  is  known. 

(b)  Available  shaft  power  is  known,  in  which  case  - 

P  *  P  ...  -  P 
c  avail  sum 


The  required  power. 


P  , 
req 


is  now  defined  as 


P  =  p  +  P 
req  sum  c 


2.7  Rate  of  Climb  Calculation 

Climb  rate  is  given  by  the  expression: 

V  -  V  siny 
v  c 


where  V  is  flight  path  velocity. 

In  axial  flight  (y  ■  0)  and,  ignoring  the  parasite  power 
component,  we  have 


2.8  Range  and  Endurance 


These  quantities  are  given  by  the  Breguet  Equations  of 
Ref.  8  namely  - 


Range 


T.V 


P  (SPC) 
sum 


1 

J 


TOR 
.. (20) 


..(21) 

.. (22) 


.. (23) 


..(24) 


.. (25) 


Equations  25  and  26  assume  that  P  /T,  V  and  SFC  axe  constants 

sum 

independent  of  the  helicopter  weight. 


3.  PROGRAM  STRUCTURE 

This  section  gives  a  description  of: 

(1)  Operating  procedure  for  running  'POLAR'. 

(2)  Basic  structure  of  calculations  made  for  a  particular 
air speed . 

(3)  Estimating  range/endurance . 

(4)  Creating  speed-power  polar. 

(5)  Estimating  power  for  a  given  V  and  y  . 

c 

(6)  Estimating  y  for  given  V  and  shaft  power. 

c 

3.1  Operating  'POLAR* 

Program  'POLAR'  provides  four  options: 

(i)  Estimate  range/endurance; 

(ii)  Provide  a  speed-power  polar  for  level  flight; 

(iii)  Estimate  the  power  required  for  a  given  airspeed  and 

rate  of  climb;  and 

(iv)  Estimate  the  possible  rate  of  climb  for  a  given 
airspeed  and  shaft  power. 

A  flow  chart  showing  the  operating  procedure  by  which  'POLAR' 
provides  these  options  is  shown  below: 


Initially  a  is  assumed  to  be  equal  to  ~Y  •  This  results  in 
quicker  convergence  of  the  iterative  calculations  and  allows  an  initial 
estimate  for  u.  X  and  6^  are  also  initially  unknown  and  are  both 
required  in  the  blade-element  relationships.  Thus  \  is  obtained  using 
simple  momentum  theory  and  is  calculated.  Using  these  values  of  y, 

X  and  9q,  a^  and  then  a  are  calculated.  This  value  of  a  is  then  used 
to  re-calculate  X  and  the  process  continues  until  X  converges. 


Having  determined  a  value  for  a,  y  is  recalculated  and 
the  process  for  determining  X,  9q,  a^^  and  a  is  repeated.  Again, 
this  loop  is  repeated  until  convergence  of  y . 

Since  values  of  the  required  quantities  have  now  been 
determined,  the  blade  angles  and  power  components  are  directly 
calculated . 


3.3  Estimating  Range  and  Endurance 


In  order  to  estimate  the  range  and  endurance ,  the  Breguet 

equations  (Eq.  25  and  26)  are  employed.  Assuming  that  SFC  is 

independent  of  velocity,  the  minimum  fuel  consumption  per  unit 

distance  and  hence  the  maximum  range  are  achieved  at  the  speed  for 

minimum  P  /V.  Similarly,  the  maximum  endurance  is  achieved  at 
sum 

the  speed  for  minimum  P 

sum 


'POLAR'  determines  range  and  endurance  in  a  manner  shown  below. 


If  range/endurance  are  required,  it  can  be  seen  from  section 

3.1  that  AV,  SFC  and  Wf  are  required.  Having  typed  in  these  parameters, 

'POLAR'  calculates  P  at  progressively  higher  velocities  until  a 

sum 

minimum  is  reached  (P  )  .  ,  from  whence  the  endurance  is  calculated. 

sum  min 

Similarly  (P  /V)  .  is  found  and  range  calculated.  Note  that  the 
sum  min 

velocity  for  maximum  endurance  is  always  lower  than  the  velocity 
for  maximum  range. 
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3.4  Creating  a  Speed-Power  Polar 

Program  'POLAR'  will  create  a  binary  data  file  which  can  be 
input  to  program  'TRANS'  (App.2)  in  order  to  obtain  graphical  or 
tabular  data  over  a  range  of  airspeeds  for  a  helicopter  in  level  flight. 
This  file  is  created  when  a  speed-power  polar  is  specified  (Section  3.1) 


The  calculations  are  accomplished  as  seen  below 


The  basic  calculations  are  made,  at  each  incremental  velocity 
specified,  starting  at  7=0  until  the  specified  7  is  reached. 

IU3JC 

3.5  Estimating  Power  for  a  Given  7elocity  and  Climb  Angle 


Power  required  for  a  given  flight  condition  is  estimated  in 
the  manner  shown  below. 


INPUT  BATA 


BASIC  CALCULATIONS 
(SECTION  3.21 


CALCULATE  Vy  <Eq.23),Pc  (Eq.20  or  24) 
Pr#q  (Eq.  22) 


OUTPUT  TO  ‘POLAR. OUT’ 


Basic  calculations  are  made  at  the  velocity  specified,  then 

V  ,  P  and  P  are  calculated  directly, 
v  c  req 

3.6  Estimating  Climb  Angle  for  Given  Velocity  and  Shaft  Power 

The  climb  angle,  y  and  thus  rate  of  climb  V  are  estimated 
c  v 

as  shown  below. 


A  first  approximation  is  made  for  y  ,  as  seen  in  Appendix  1, 

and  the  basic  calculations  are  made  using  this  value  for  y  .  Pc  and 

then  y  are  calculated.  If  this  value  of  y  does  not  agree  with  the 
c  c 

first  approximation,  then  this  value  is  used  in  the  basic  calculations 
and  the  process  repeated  until  converges. 

4.  EXAMPLE  OP  PERFORMANCE  ESTIMATION 


1. 


Speed-Power  Polar  in  level  flight. 


2.  Range  and  Endurance. 

3 .  Climb  Performance . 

3.1  Rate  of  climb  possible  for  a  given  airspeed 
and  shaft  power. 

3.2  Shaft  power  required  for  a  given  airspeed  and 
rate  of  climb. 

3.3  Maximum  rate  of  climb. 

3.4  Power  required  for  a  given  vertical  climb. 

4.  Ceiling  Estimates. 

4.1  Absolute  ceiling. 

4.2  Service  ceiling. 

4.3  Hover  ceiling  out  of  ground  effect. 

4.4  Hover  ceiling  in  ground  effect. 

The  helicopter  is  assumed  to  have  the  following 
specifications : 

W  «  9604  lb  (42718  N) 

Wf  »  1439  lb  (6400  N) 

OR  »  717.2  ft/s  (218.6  m/s) 

R  «  21  ft  (6.4  m) 

c  »  1.3  ft  (0.395  m) 

b  -  4 

6  *  -8° 

Y  =  15 

B  -  0.97 

a  -  5.73/rad 

f  =  17  fta  (1.58  ma) 


P. 

inst 


1500  Hp 


SFC-  1.11  lb/Hp-Hr  (0.504  kg/Hp-Hr) 


The  main  rotor  blades  fore  assumed  to  be  of  NACA  23012 


aerofoil  section,  with  drag  coefficient  related  to  angle  of  attack  by 

the  following  expression:  *  0.0087  -  0.0216  ar  +  0.4  a  2 

o 

therefore  6^  *  0.0087,  6.,  ■  -0.0216,  6^  “  0.4. 

It  is  also  assumed  that  the  helicopter  is  operating  at 
ISA  +  15°C  conditions,  unless  otherwise  stated.  (Appendix  3) 

Appendix  1  shows  how  program  'POLAR'  is  loaded  onto  a  computer 
and  explains  how  running  ' POLAR '  creates  either  a  ' . DAT '  file  or  a 
’.OUT*  file.  A  '.DAT'  file  requires  further  processing  by  program 
'TRANS'  (Appendix  2)  to  give  column  and  graphical  output. 

4.1  Speed-Power  Polar  in  Level  Flight 

For  the  example  following,  metric  units  cure  incorporated 
for  the  helicopter  input  data.  When  a  Speed-Power  polar  is  specified, 
the  results  are  obtained  over  a  range  of  speeds  and  output  in  binary 
form  to  file  'POLAR.DAT' .  'TRANS'  can  then  be  run  to  obtain  a  plot 
of  power  against  speed  (Fig.  3).  Alternatively,  'TRANS'  may  be  run 
to  produce  a  column  print  of  the  same  information,  as  shown  in  Appendix  2 
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.RU  POLAR 

TITLE  (TWO  LINES  OF  UP  TO  48  CHARACTERS) 

:  BLADE  ELEMENT  METHOD 

POWER  SPEED  POLAR 

ARE  UNITS  IN  IMPERIAL  OR  METRIC  ( I  OR  M)  ?  M 

ALL  UP  WEIGHT  CN  OR  LB)  -  ?  42718 

ROTOR  TIP  SPEED  CM/S  OR  FT/S)  -  ?  218.6 

ROTOR  RADIUS  CM  OR  FT)  -  ?  6.4 

ROTOR  BLADE  CHORD  CM  OR  FT)  -  ?  0.395 

NUMBER  OF  ROTOR  BLADES  -  ?  4 

2D  LIFT  CURVE  SLOPE  (/RAO)  -  ?  5.73 

TIP  LOSS  FACTOR  -  ?  0.97 

BLADE  TWIST  (DEG)  -  ?  -8 

EQUIVALENT  FLAT  PLATE  AREA  (M**2  OR  FT**2 )  -  ?  1.58 
LOCKS  NUMBER  -  ?  15 

DRAG  POLAR  CO-EFFICIENT  (DELO)  -  ?  0.0087 
OELTA1  (/RAO)  -  ?  -0.0216 
DELTAS  ( /PAD-* ■*'2 )  -  ?  0.4 

AIR  DENSITY  (KG/M**3  OR  SLUG/FT*»3 )  *  ?  1 . 16** 

IS  RANGE/ENDURANCE  REQUIRED  ?  (Y  OR  N)  :  N 
IS  A  SPEED-POWER  POLAR  REQUIRED  ?  (Y  OR  N)  :  Y 
MAXIMUM  VELOCITY  (KNOTS)  -  ?  150 
VELOCITY  INTERVAL  (KNOTS)  -  ?  2 

END  OF  EXECUTION 

CPU  TIME:  6.10  ELAPSED  TIME:  53.40 
EXIT 


4.2  Range  and  Endurance 

The  example  below  shows  how  'POLAR'  is  run  to  obtain  the 
range  of  the  hypothetical  helicopter.  The  data  is  input  in  metric 
units  and  the  results  are  output  to  file  'POLAR. OUT'.  'POLAR'  is 
run  as  in  the  previous  example,  with  the  following  changes. 
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IS  RANGE/ENDURANCE  REQUIRED  ?  (Y  OR  N)  :  Y 
UELOCITY  INTERVAL  (KNOTS)  -  ?  2.0 

SPECIFIC  FUEL  CONSUMPTION  ?  (KG/HP- HR  OR  LB/HP- HR) 
UEIGHT  OF  FUEL  CARRIED  ?  (N  OR  LB)  :  6400 

STOP 

ENO  OF  EXECUTION 

CPU^TIME:  3.60  ELAPSED  TIME:  55.02 


.TY  POLAR. OUT 
&LADE  ELEMENT  METHOD 
RANGE/ENDURANCE  ESTIMATE 
ALL  UP  UEIGHT  -  42718.0  N 

AIR  DENSITY  -  1.1644  KG/M**3 

TIP  SPEED  -  218.6  M/'S 

ROTOR  RADIUS  -  6.4  M 

ROTOR  CHORD  -  0.395  M 

NUMBER  OF  ROTOR  BLAOES 
LIFT  CURVE  SLOPE 
TIP  LOSS  FACTOR  - 
BLADE  TUI  ST  -  -8.0  DEG 

EQUIVALENT  FLAT  PLATE  AREA  »  1.58  M**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT(DELTA  0)  -  0.0087 

OELTA  1  -  -0.0216  /RAD 

DELTA  2  -  0.4000  /RAD**2 

S.F.C.  -  0.504  KG/HP-HR 

FUEL  UEIGHT  -  6400.0  N 


0.504 


4.0 
5.73  /RAO 
0.97 


RESULTS  FOR  ENDURANCE  OF  2.81  HOURS,  ARE  AS  FOLLGUS  : 

AIRSPEED  -  72.00  KNOTS 

ADVANCE  RATIQ(MU)  -  0.169 

I NOUCED  VELQC I TY ( NU )  -  3.30  M/S 

INFLOU  RAT  I G( LAMBDA)  -  -0.03  06 

FLAT  PLATE  ORAG  -  1262.0  N 

THRUST  -  42736.6  N 

THRUST  CO-EFFICIENT  <CT)  -  0.00597 

BLAOE  SOLIDITY  -  0.0786 

COLLECTIVE  (THETA  0)  -  13.1  DEG 

CONING  ANGLE  (aO)  -  7.8  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (al)  -  2.8  DEG 

LATERAL  FLAPPING  ANGLE  (bl)  -  1.8  DEG 

OI SC  ANGLE  OF  ATTACK  -  -4.5  DEG 

INDUCED  POUER  -  217.8  HP 

PARASITE  POUER  -  62.7  HP 

PROFILE  POUER  -  220.1  HP 

TOTAL  POUER  «  500.6  HP 

LOCAL  ANGLE  OF  ATTACK ( 027ODEG&4OS  RADIUS  )  -  7.2  DEG 

LOCAL  ANGLE  OF  ATTACK ( 027ODEG&T I P )  -  6.4  DEG 

LOCAL  ANGLE  OF  ATTACK  ( 027QDEG&PT . OF  UT-0.4)  -  8.2  DEG 

RADIAL  POINT  OF  UT-0.4  (@270DEG)  -  0.57 


1  RESULTS  FOR  RANGE  OF 


471.5  KM,  ARE  AS  FOLLOUS 


AIRSPEED  -  112.00  KNOTS 

ADVANCE  RATIO (MU)  -  0.261 

INDUCED  VELOCITY(NU)  -  2.46  M/S 

INFLOU  RAT IO (LAMBDA )  -  -0.0510 

FLAT  PLATE  DRAG  -  3053.8  N 

THRUST  -  42827.0  N 

THRUST  CO-EFFICIENT  (CT)  -  0.00598 

BLADE  SOLIDITY  -  0.0786 

COLLECTIVE  (THETA  0)  -  14.3  DEG 

CONING  ANGLE  CaO)  -  7.8  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (al) 


4.6  DEG 


LATERAL  FLAPPING  ANGLE  (bl)  -  2.7  DEG 

DISC  ANGLE  OF  ATTACK  -  -3.7  DEG 

INDUCED  POUER  -  141.4  HP 

PARASITE  POUER  -  236.0  HP 

PROFILE  POUER  -  241.7  HP 

TOTAL  POUER  -  619.1  HP 

LOCAL  ANGLE  OF  ATTACK ( @2 700EG&40K  RADIUS  )  - 
LOCAL  ANGLE  OF  ATTACK (B270DEG4T IP)  -  8.5  DEG 

LOCAL  ANGLE  OF  ATTACK  (027QOEG4PT. OF  UT-0.4)  - 


3.3  DEG 


9.2  DEG 


4.3  Climb  Performance 


In  order  to  estimate  climb  performance  using  'POLAR' ,  the 
power  available  at  the  main-rotor  shaft  must  be  known.  The  installed 
engine  power  of  the  hypothetical  helicopter  is  given  as  1500  Hp. 

However,  not  all  of  this  power  is  available  at  the  main  rotor  shaft. 

Some  power  will  be  absorbed  by  the  main-rotor  transmission,  accessories 
(hydraulics,  pneumatics,  electronics,  heating  and  cooling  systems  etc.) 
and  the  tail-rotor.  We  will  assume  that  the  sum  of  these  power  losses 
is  15%  of  the  installed  engine  power. 

It  should  be  noted  that  in  low  speed  forward  flight,  the  power 
absorbed  by  the  tail-rotor  is  proportional  to  the  power  required  by  the  main 
rotor  which  in  turn  varies  with  forward  speed  (Fig.  3) .  Thus  in  reality 
the  power  losses  will  not  be  constant  at  various  velocities.  (Ref.  6 
states  that  tail-rotor  losses  vary  from  6%  of  total  power  in  hover, 
to  3%  in  forward  flight) .  For  simplicity  we  will  assume  15%  losses 
throughout  the  speed  range.  Therefore,  the  available  power  at  the  main 
rotor  shaft  is  1275  Hp. 

It  should  also  be  noted  that  in  fast  descents  (negative  climbs) . 
where  |Vv|  i  v,  helicopters  may  enter  the  'Vortex  Ring  State'.  When 
this  phenomenon  occurs,  momentum  theory  breaks  down  and  'POLAR'  will 
not  give  meaningful  results. 

4.3.1  Rate  of  climb  for  a  given  airspeed  and  shaft  power 


In  the  example  below,  'POLAR'  is  used  to  estimate  the  rate 
of  climb  at  a  flight  path  airspeed  of  50  kt.  The  input  data  is  in 
Imperial  units. 


RU  POLAR 


TITLE  (TWO  LINES  OF  UP  TO  60  CHARACTERS) 

:  BLADE  ELEMENT  METHOD  -  IMPERIAL  UNITS 

RATE  OF  CLIMB  REQUIRED  FOR  GIUEN  AIRSPEED  &  SHAFT  POWER 

ARE  UNITS  IN  IMPERIAL  OR  METRIC  (I  OR  M)  ?  I 

ALL  UP  WEIGHT  (N  OR  LB)  -  ?  9604 

ROTOR  TIP  SPEED  (M/S  OR  FT/S)  -  ?  717.2 

ROTOR  RADIUS  (M  OR  FT)  -  ?  21 

ROTOR  BLADE  CHORD  (M  OR  FT)  -  ?  1.3 

NUMBER  OF  ROTOR  BLADES  -  ?  4 

2D  LIFT  CURUE  SLOPE  (/RAO)  -  ?  5.73 

TIP  LOSS  FACTOR  -  7  0.97 

SLADE  TWIST  (DEG)  -  ?  --3 

EQUIVALENT  FLAT  plats  AREA  ( CR  F 7**2 )  -  ?  17 
LOCKS  NUMBER  -  ?  15 

DRAG  POLAR  CO-EFFICIENT  (DELO)  -  ?  0.0097 
OELTA1  (/RAO)  -  ?  -0.0216 
DELTA2  (/RAO*+2)  -  ?  0.4 

AIR  DENSITY  (KG/M**3  OR  SLUG/FT**3 )  -  ?  0.002259 

IS  RANGE/ENDURANCE  REQUIRED  ?  (Y  OR  N)  :  N 

IS  A  SPEED- POWER  POLAR  REQUIRED  ?  (Y  OR  N)  :  N 

IS  POWER  REQUIRED  FOR  A  GIUEN  CLIMB  ANGLE 
i.e.  RATE  OF  CLIMB  (Y  OR  N)  ?  N 

UELOCITY  IN  KNOTS  ?  :  50. 

WHAT  IS  POWER  AVAILABLE  AT  MAIN  ROTOR  SHAFT  (HP)  ?  :  1275 
STOP 

END  OF  EXECUTION 

CPU  TIME:  1.65  ELAPSED  TIME:  1:6.02 
EXIT 
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The  results  of  this  run,  are  output  onto  File  'POLAR. OUT' 


.TY  POLAR. OUT 

&LAOE  ELEMENT  METHOD  -  IMPERIAL  UNITS 

RATE  OF  CLIMB  REQUIRED  FOR  GIUEN  AIRSPEED  &  SHAFT  POWER 
ALL  UP  WEIGHT  -  9604.0  LB 

AIR  DENSITY  -  0.002259  SLUG/FT**3 

TIP  SPEED  -  717.2  FT/'S 

ROTOR  RADIUS  -  21.0  FT 

ROTOR  CHORD  -  1 . 3  FT 

NUMBER  OF  ROTOR  SLADES  -  4.0 

LIFT  CURUE  SLOPE  -  5.73  /RAD 

TIP  LOSS  FACTOR  -  0.97 

BLADE  TWIST  -  -8.0  DEG 

EQUIUALENT  FLAT  PLATE  AREA  -  17.0  FT**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

DELTA  1  -  -0.0216  /RAO 
DELTA  2  -  0.4000  /RAO**2 

AIRSPEED  -  50.00  KNOTS 


0  .  1 0  0 

-  17.33  FT/S 

-  -0.0875 

7  !  3 


.JO 


ADUANCE  RAT  10 (MU)  - 
INDUCED  UELOCITY(NU) 

INFLOW  RAT IO( LAMBDA) 

-LA 7  PLATE  DRAG  - 
T'-RLST  *  ro7"..  1  _3 

THRUST  CO- EFF I C I  ENT  C C7 )  »  3 .  3  0  o  3 1 

SLAOE  5CL 1 0  I TY  -  0.0738 

COLLECT  I UE  K THETA  0)  -  13.3  OEG 

CONING  ANGLE  (aO)  *  3.9  OEG 

LONGITUDINAL  FLAPPING  ANGLE  Cal)  -  2.4  DEG 

LATERAL  FLAPPING  ANGLE  <bl>  -  1.2  DEG 

DISC  ANGLE  OF  ATTACK  -  -32.1  DEG 
INDUCED  POWER  -  313.5  HP 

PARASITE  POWER  -  21.0  HP 

PROFILE  POWER  -  225.1  HP 

TOTAL  POWER  -  559.6  HP 

REQU I RED/AU A I LABLE  SHAFT  POWER  -  1275.0  HP 

CLIMB  POWER  -  715.4  HP 

RATE  OF  CLIMB  -  2458.1  FT/M  I N 

CLIMB  ANGLE  -  29.0  DEG 

LOCAL  ANGLE  OF  ATTACK (02 7 0DEG2,40J<  RADIUS  )  - 
LOCAL  ANGLE  OF  ATTACK ( B2 70 DEG&T IP)  -  7.4  DEG 

LOCAL  ANGLE  OF  ATTACK  I 027ODEG&PT . OF  UT-0.4)  - 
RAO  I AL  POINT  OF  UT-0.4  C@270DEG)  -  0.500 


1.8  DEG 
4.8  DEG 


From  the  previous  example  we  can  see  that  at  an  airspeed  of 
50  kts,  power  of  1275  Hp,  the  rate  of  climb  is  2458.1  ft/min  and  climb 
angle  is  29°. 

The  inverse  case  is  shown  in  the  example  below,  where  ‘POLAR 
is  used  to  estimate  the  power  required  at  50  kt  and  climb  angle  of  29® 
'POLAR*  is  run  as  in  the  previous  example,  with  the  following  changes. 

IS  POWER  REQUIRED  FOR  A  GIVEN  CLIMB  ANGLE 

i.e.  RATE  OF  CLIMB  (Y  OR  N)  ?  Y 

VELOCITY  IN  KNOTS  ?  :  50. 

WHAT  IS  CLIMB  ANGLE  -  GAMMA  -  ASINCUU/U)  -  (DEG)  ?  :  29 

STOP 

END  OF  EXECUTION 

CPU  TIME:  1.44  ELAPSED  TIME:  50.60 

EXIT 


Results  typed  below  agree  with  those  of  the  previous  example 
Any  small  differences  result  from  rounding  errors. 


.TY  POLAR. OUT 

BLADE  ELEMENT  METHOD  -  IMPERIAL  UNITS 
POWER  REQUIRED  FOR  A  GIUEN  AIRSPEED  &  CLIMB  ANGLE 
ALL  UP  WEIGHT  -  9604.0  LB 

AIR  DENSITY  -  0.002259  SLUG/FT**3 

TIP  SPEED  -  717.2  FT/S 

ROTOR  RADIUS  -  21.0  FT 

ROTOR  CHORD  -  1 . 3  FT 

NUMBER  OF  ROTOR  BLADES  -  4.0 

LIFT  CURUE  SLOPE  -  5.73  /RAD 

TIP  LOSS  FACTOR  «  0.97 

BLADE  TWIST  -  -9.0  DEG 

EQUIUALENT  FLAT  PLATE  AREA  -  17.0  FT**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

DELTA  1  -  -0.0216  /RAD 
DELTA  2  -  0.4000  /RAO*+2 

AIRSPEED  -  50.00  KNOTS 

AOUANCE  RATIO(MU)  -  0.100 

INDUCED  UELOC ITY(NU)  -  17.33  FT/S 

INFLOW  RATIO(LAMBOA)  -  -0.097* 

FLAT  PLATE  DRAG  -  136.7  L3 

THRUST  -  3-71.0  •_  = 

THRUST  CC-Sr “ ' C IENT  tCT)  -  0.00601 

3LADE  SOLIDITY  -  0.0793 

COLLECT  I UE  (THETA  0)  -  13.3  OEG 

CONING  ANGLE  (aO)  -  9.9  DEG 

LONGITUDINAL  FLAPPING  ANGLE  Cal)  -  2.4  DEG 

LATERAL  FLAPPING  ANGLE  <bl)  -  1.2  DEG 

DISC  ANGLE  OF  ATTACK  -  -32.1  DEG 

INDUCED  POWER  -  313.5  HP 

PARASITE  POWER  -  21.0  HP 

PROFILE  POWER  -  225.1  HP 

TOTAL  POWER  -  559.6  HP 

REQU I RED/AUA I LABLE  SHAFT  POWER  -  1274.0  HP 

CLIMB  POWER  -  714.4  HP 

RATE  OF  CLIMB  -  2454.8  FT/MIN 

CLIMB  ANGLE  -  29.0  DEG 

LOCAL  ANGLE  OF  ATTACK ( 027QDEG8.4OS  RADIUS  )  -  1.9  DEG 

LOCAL  ANGLE  OF  ATTACK ( 027ODEG&T I P )  -  7.4  DEG 

LOCAL  ANGLE  OF  ATTACK  I  (3270DEG&PT . OF  UT-0.4)  -  4.8  DEG 

RADIAL  POINT  OF  UT-0.4  C027ODEG)  -  0.500 
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4.3.3  Maximum  rate  of  climb 

The  maximum  rate  of  climb  occurs  when  the  power  required  for 
level  flight  is  a  minimum.  Using  results  of  §4.2  or  figure  3,  72  kt 
is  the  velocity  at  which  maximum  climb  rate  will  be  achieved.  This  is 
demonstrated  in  the  example  below.  'POLAR*  is  run  as  in  §4.3.1  but 
with  a  velocity  of  72  knots.  Only  the  output  file  is  shown. 


.TY  POLAR. OUT 

Slade  element  method  -  imperial  units 

MAXIMUM  RATE  OF  CLIMB 
ALL  UP  WEIGHT  -  9604.0  LB 

AIR  DENSITY  -  0.002259  SLUG/FT**? 

TIP  SPEED  -  717.2  FT/S 

ROTOR  RAD  I  US  -  21.0  FT 

ROTOR  CHORD  -  1 . 3  FT 

NUMBER  OF  ROTOR  BLA0E5  -  4.0 

LIFT  CUPUE  SLOPE  -  5.73  /RAD 

TIP  LOSS  FACTOR  -  0.97 

BLADE  TWIST  -  -3.0  DEG 

EQUIVALENT  "LA7  PLATE  AREA  -  17.0  FT*-* 2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CC-SFF; O'ENT  » DELTA  3;  -  0.0037 

DELTA  1  -  -0.0216  /RAO 
DELTA  2  -  0.4000  /RAO**2 

AIRSPEED  -  72.00  KNOTS 

ADVANCE  RAT  1 0 ( MU )  -  0.152 

INDUCED  VELOCITYCNU)  -  12.13  FT/S 

INFLOW  RAT  1 0 ( LAMBDA )  -  -0.0923 

FLAT  PLATE  DRAG  -  233.6  LB 

THRUST  -  9710.1  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.00603 

BLAOE  SOLIDITY  -  0.0788 

COLLECTIVE  (THETA  0)  *  18.4  DEG 

CONING  ANGLE  (aO)  -  3.8  DEG 

LONG  I TUO INAL  FLAPPING  ANGLE  (al)  -  3.7  DEG 

LATERAL  FLAPPING  ANGLE  (bl)  -  1.8  DEG 

DISC  ANGLE  OF  ATTACK  -  -26.4  DEG 

INDUCED  POWER  -  215.0  HP 

PARASITE  POWER  -  62.7  HP 

PROFILE  POWER  -  231.9  HP 

TOTAL  POWER  -  509.5  HP 

REQUIRED/AVAILABLE  SHAFT  POWER  -  1275.0  HP 

CLIMB  POWER  -  765.5  HP 

RATE  OF  CLIMB  -  2630.2  FT/MIN 

CLIMB  ANGLE  -  21.1  DEG 

LOCAL  ANGLE  OF  ATTACK(@270DEG&40S  RADIUS  )  - 
LOCAL  ANGLE  OF  ATTACK ( §270 DEGiT I P )  -  8.5  DEG 

LOCAL  ANGLE  OF  ATTACK  (027ODEG&PT , OF  UT-0.4)  - 
RADIAL  POINT  OF  UT-0.4  (@270DEG)  -  0.552 


0.2  DEG 
5.9  DEG 


Vertical  climb  performance  can  be  estimated  by  setting  the 
climb  angle  to  90°  and  specifying  the  rate  of  climb.  If  the  vertical 
rate  of  climb  is  required  for  a  given  power,  the  program  must  be  run 
several  times.  Each  time  the  flight  path  velocity  must  be  varied  until 
the  power  is  equal  to  that  available  in  the  helicopter. 


The  example  below  shows  the  final  result  for  the  maximum 
vertical  rate  of  climb. 


§LADE°kLEM?NT  METHOD  -  IMPERIAL  UNITS 
POWER  REQUIRED  FOR  A  UERTICAL  CLIMB 
ALL  UP  WEIGHT  -  9604.0  LB 

AIR  DENSITY  -  0.002259  SLUG/FT**3 

TIP  SPEED  -  717.2  FT/S 

ROTOR  RADIUS  -  21.0  FT 

ROTOR  CHORD  -  1 . 3  FT 

NUMBER  OF  ROTOR  BLADES  -  4.0 

LIFT  CURUE  SLOPE  -  5.73  /RAD 

TIP  LOSS  FACTOR  -  0.97 

BLADE  TWIST  -  -0.0  DEG 

EQUIUALENT  FLAT  PLATE  AREA  -  17.0  FT**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

DELTA  1  -  -0.0216  /RAO 
DELTA  2  -  0.4000  /RAO**2 

AIRSPEED  -  11.30  KNOTS 

AOUANCE  RATIO (MU)  -  0.000 

INDUCED  UELQCITY(NU)  -  30.79  FT/S 

INFLOW  RAT  1 0 ( LAMBDA )  -  -0.0842 

FLAT  PLATE  ORAG  -  9 . 0  LB 

THRUST  -  9613.0  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.00597 

BLADE  SOLIDITY  -  0.0788 

COLLECT  I UE  (THETA  0)  -  18.2  DEG 

CONING  ANGLE  (aO)  -  8.9  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (al)  -  0.0  DEG 

LATERAL  FLAPPING  ANGLE  (bl)  -  0.0  DEG 

DISC  ANGLE  OF  ATTACK  -  -90.0  DEG 

INDUCED  POWER  -  678.0  HP 

PARASITE  POWER  -  0 . 4  HP 

PROFILE  POWER  -  219.3  HP 

TOTAL  POWER  -  897.6  HP 

REQU I RED/AUA I LABLE  SHAFT  POWER  -  1274.8  HP 

CLIMB  POWER  -  377.2  HP 

RATE  OF  CLIMB  -  1296.2  FT/M  IN 

CLIMB  ANGLE  -  90.0  DEG 

LOCAL  ANGLE  OF  ATTACK ( 027ODEG&4OK  RADIUS  )  -  3.1  DEG 

LOCAL  ANGLE  OF  ATTACK(@270DEG2.TIP)  -  5.4  DEG 

LOCAL  ANGLE  OF  ATTACK  (@270DEG&PT. OF  UT-0.4)  -  3.1  DEG 

RADIAL  POINT  OF  UT-0.4  (@270DEG)  -  0.400 


As  can  be  seen  above,  the  power  required  for  a  vertical  climb 
rate  of  1296.2  ft/min  (1274. S  Hp)  is  almost  equal  to  the  power 
available  (1275  Hp) . 
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4.4  Ceiling  Estimates 

The  helicopter  performance  analyst  is  interested  in  four 
ceilings ,  namely : 

(1)  Absolute  Ceiling  -  defined  as  the  maximum  altitude  at 
which  the  climb  rate  is  zero. 

(2)  Service  Ceiling  -  defined  as  the  maximum  altitude  at 
which  the  climb  rate  is  100  ft/min  (0.5  m/s). 

(3)  Hover  Ceiling  O.G.E.  -  defined  as  the  maximum  altitude 
at  which  the  power  available  equals  the  power  required 
to  hover  out  of  ground  effect. 

(4)  Hover  Ceiling  I.G.E.  -  defined  as  the  maximum  altitude 
at  which  the  power  available  equals  the  power  required 
to  hover  in  ground  effect. 

Reference  9  presents  curves  of  engine  power  variation  with 
altitude,  which  can  be  approximated  by  the  following  relationship: 


P 


avail 


P„  (1.13  A  -  0.13) 
SL 


where  P  is  the  shaft  power  available  at  sea-level  (1275  Hp) 
SL 


and 


A  =  p/p 


SLISA+15 


[It  should  be  noted  that  actual  engine  power  variation  with  altitude 
should  be  used,  if  known] . 

Variations  of  density  with  atmospheric  conditions  may  be 
found  using  the  locally  available  program  'ATMOS'  (see  Appendix  3). 

At  high  altitudes,  conditions  are  such  that  the  main  rotor 

may  operate  at  blade  angles  of  attack  where  partial  stall  may  occur. 

Gustafson  and  Gessow  (Ref.  10)  found  by  flight  tests  that  the  ratio  of 

measured  profile  power  (P  )  to  calculated  profile  power  (P  )  was 

meas  calc 

unity  until  an  „  reached  a  value  of  12°.  This  ratio  rose  linearly 

(J,  /UJ 


reached 


with  aM  to  a  value  of  approximately  two  when  a,, 

A/UJ  (1 1  2.  IQ) 

16°.  (see  diagram  below) .  Above  this  value  the  effects  of  stall 

were  such  that  the  helicopter  became  too  difficult  to  control. 


P„ 

°ni*as 


,270) 


4.4.1  Absolute  Ceiling 

Absolute  ceiling  will  be  attained  when  the  minimum  cower  in 

level  flight  is  equal  to  the  power  available.  'POLAR'  is  ran  at  a 

series  of  altitudes  (densities)  using  the  RANGE /ENDURANCE  option, 

since  minimum  power  is  calculated  in  order  to  estimate  endurance. 

[Note:  If  and  S.F.C.  are  unknown,  arbitrary  values  may  be  assigned 

to  these  parameters,  since  they  are  only  used  to  calculate  range/ 

endurance  and  have  no  bearing  on  the  power  calculations] .  At  each 

altitude  the  minimum  power  is  obtained.  The  value  of  210) 

then  examined  and,  if  necessary,  ?  is  adjusted  to  allow  for  stall  as 

o 

stated  previously.  Thus  the  true  total  power  is  obtained  and  then 
compared  with  the  power  available  at  that  altitude.  The  process 
described  above  is  repeated  until  convergence  is  achieved. 

The  example  below  is  the  final  result  obtained  after  several 
iterations  and  corresponds  to  an  altitude  of  17300  ft  at  ISA  +  15 °C 
conditions . 
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.RU  POLAR 


TITLE  (TWO  LINES  OF  UP  TO  60  CHARACTERS) 

:  SLADE  ELEMENT  METHOD 

ABSOLUTE  CEILING  ESTIMATION 

ARE  UNITS  IN  IMPERIAL  OR  METRIC  (I  OR  M)  ?  I 

ALL  UP  WEIGHT  (N  OR  LB)  -  ?  9604 

ROTOR  TIP  SPEED  (M/S  OR  FT/S)  -  ?  717.2 

ROTOR  RADIUS  CM  OR  FT)  -  ?  21 

ROTOR  BLADE  CHORD  (M  OR  FT)  -  ?  1.3 

NUMBER  OF  ROTOR  BLADES  -  ?  4 

2D  LIFT  CURUE  SLOPE  (/RAD)  -  ?  5.73 

TIP  LOSS  FACTOR  -  ?  0.97 

BLAOE  TWIST  (DEG)  -  ?  -8 

EQUIUALENT  FLAT  PLATE  AREA  (M**2  OR  FT**2 )  -  ?  17 
LOCKS  NUMBER  -  ?  15 

DRAG  POLAR  CO-EFFICIENT  (OSLO)  -  ?  0.0087 
DELTA1  (/RAO)  -  ?  -0.0216 
DELTA2  (/RAO**2)  -  ?  0.4 

AIR  DENSITY  (KG/M**3  OR  SLUG/FT**3 )  -  ?  0.001309 
IS  RANGE/ENDURANCE  REQUIRED  ?  (Y  OR  N)  :  Y 
UELOCITY  INTERUAL  (KNOTS)  -  ?  2 

SPECIFIC  FUEL  CONSUMPTION  ?  (KG/HP- HR  OR  LB/HP-HR)  :  1.11 
WEIGHT  OF  FUEL  CARRIED  ?  (N  OR  LB)  :  1439 

STOP 

END  OF  EXECUTION 

CPU  TIME:  4.99  ELAPSED  TIME:  1:56.78 
EXIT 


TV  POLAR. OUT 

BLADE  ELEMENT  METHOD 
ABSOLUTE  CEILING  ESTIMATION 
ALL  UP  WEIGHT  -  9604.0  LB 

AIR  DENSITY  -  0.001309  SLUG/FT**? 

TIP  SPEED  -  717.2  FT/S 

ROTOR  RADIUS  -  21.0  FT 

ROTOR  CHORD  -  1 . 3  FT 

NUMBER  OF  ROTOR  BLADES  -  4.0 

LIFT  CURUE  SLOPE  -  5.73  /RAD 

TIP  LOSS  FACTOR  -  0.97 

BLADE  TWIST  -  -0.0  DEG 

EQU I UALENT  FLAT  PLATE  AREA  -  17.0  FT**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

DELTA  1  -  -0.0216  /RAO 

DELTA  2  -  0.4000  /RAD**2 

S.F.C.  -  1.11  LB/HP-HR 

FUEL  WEIGHT  -  1439.0  LB 


RESULTS  FOR  ENDURANCE  OF  2.49  HOURS,  ARE  AS  FOLLOWS  : 

AIRSPEED  -  94.00  KNOTS 

ADUANCE  RAT IO (MU)  -  0.219 

INDUCED  UELOCITY(NU)  -  16.38  FT/S 

INFLOW  RAT  I O ( LAMBDA )  -  -0.0530 

FLAT  PLATE  DRAG  -  280.1  LB 

THRUST  -  9608.1  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.01030 

BLAOE  SOLIDITY  -  0.0788 

COLLECT  I UE  (THETA  0)  -  13.1  DEG 

CCNISiG  ANGLE  (a0)  -  13.3  DEG 

-ENG  :  ~-0  :  NAL  ~L.-F°  :  MG  -NGLE  ( a  1 )  -  3  .  2  "EG 

LATERAL  FLAP® ING  ANGLE  CbiJ  -  4.0  DEG 

DISC  ANGLE  CF  ATTACK  -  -7.3  DEG 

INDUCED  POWER  -  286.2  HP 

PARASITE  POWER  -  30.8  HP 

PROFILE  POWER  -  196.7  HP 

TOTAL  POWER  -  563.7  HP 

LOCAL  ANGLE  OF  ATTACK ( (327 0DEG&40S  RADIUS  )  -  11.7  DEG 


__  _ _  .  11.7  DEG 

LOCAL  ANGLE  OF  ATTACK ( 3270 DEG&T I P )  -  14.1  DEG 

LOCAL  ANGLE  OF  ATTACK  T3270DEG&PT. OF  UT-0.4)  -  15.1  DEG 

RAOIAL  POINT  OF  UT-0.4  (@270OEG)  -  0.619 

1  RESULTS  FOR  RANGE  OF  335.4  MILES,  ARE  AS  FOLLOWS  : 

AIRSPEED  -  142.00  KNOTS 

ADUANCE  RATIO (MU)  -  0.325 

INDUCED  UELOCITY(NU)  -  10.94  FT/S 

INFLOW  RAT  I Q( LAMBDA)  -  -0.0935 

FLAT  PLATE  DRAG  -  639.1  LB 

THRUST  -  9625.2  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.01032 

BLADE  SOLIDITY  -  0.0788 

COLLECT I UE  (THETA  0)  -  20.1  DEG 

CONING  ANGLE  (aO)  «  13.6  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (al)  -  9.7  DEG 

LATERAL  FLAPPING  ANGLE  (bl)  -  5.8  DEG 

DISC  ANGLE  OF  ATTACK  -  -13.5  DEG 

INDUCED  POWER  -  191.5  HP 

PARASITE  POWER  -  278.5  HP 

PROFILE  POWER  -  215.7  HP 

TOTAL  POWER  -  685.7  HP 

LOCAL  ANGLE  OF  ATTACK ( 3270DEG&40S  RADIUS  )  -  -1.9  DEG 

LOCAL  ANGLE  OF  ATTACK ( S270DEG&T I P )  -  18.5  DEG 

LOCAL  ANGLE  OF  ATTACK  (3270DEG&PT. OF  UT-0.4)  «  18.5  DEG 

RAOIAL  POINT  OF  UT-0.4  (@270DEG)  -  0.725 
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From  these  results,  the  minimum  power  cccurs  at  94  knots: 

f total  ’  563'7  «■'  Po  ,  '  196-7  Hp'  “(1,270)  ■  14‘1* 

calc 

Now, 


and 


req 


total 


+  P 


meas 


P 

°calc 


thus 


P  -  P  +  (  - —  )  P 

o  o  ,  x  4°  o 

meas  calc  calc 


p  ,  p  +  (  a(l,270)  12*}  p 

req  total  v  4°  ’  o 


calc 


563.7  +(  14..-.1-f.1.2°)  196.7 


*  667  Hp 


Power  available  at  the  shaft  is  given  by 

P  .  =  1275  (1.13x0.001309/0.002259-0.13) 

avail 

=  669  Hp 


Since  the  power  available  is  approximately  equal  to  the  power 
required,  the  absolute  ceiling  is  approximately  17,300  ft. 


4.4.2  Service  ceiling 

The  climb  rate  for  this  case  is  defined  to  be  100  ft/min 
(1.67  ft/s).  As  for  4.4.1,  ’POLAR'  must  be  run  for  a  series  of  altitudes 
(densities)  for  level  flight  conditions.  Having  determined  the  velocity 
for  minimum  total  power,  y can  be  calculated  from  equation  (23).  'POLAR' 
cam  then  be  run  using  the  specified  yc  and  V.  This  process  is  continued 
until  convergence  of  total  power  required  and  power  available. 


The  example  below  shows  how  to  estimate  service  ceiling  and 
is  the  final  result  after  several  iterations. 

'POLAR'  is  run  as  for  the  previous  example,  to  find  the  velocity 
for  minimum  power. 

The  output  file  shown  below  corresponds  to  an  altitude  of  16900  ft 


,  TY  POLAR. OUT 

SERVICE  CEILING 
VELOCITY  £  MINIMUM  POWER 
ALL  UP  WEIGHT  -  9604.0  LB 

AIR  OENSITY  -  0.001327  SLUG/FT**3 

TIP  SPEED  -  717.2  FT/S 

ROTOR  RADIUS  -  21.0  FT 

ROTOR  CHORD  -  1 . 3  FT 

NUMBER  OF  ROTOR  BLADES  -  4.0 

LIFT  CURVE  SLOPE  -  5.73  /RAD 

TIP  LOSS  FACTOR  -  0.97 

BLADE  TWIST  -  -0.0  DEG 

EQUIVALENT  FLAT  PLATE  AREA  -  17.0  FT**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

DELTA  1  -  -0.0216  /RAD 

OELTA  2  -  0.4000  /RAD**2 

S.F.C.  -  1.11  LB/HP-HR 

FUEL  WEIGHT  -  1439.0  LB 

RESULTS  FOR  ENDURANCE  OF  2.50  HOURS,  APE  AS  FOLLOWS  : 

AIRSPEED  -  92.00  KNOTS 

AOVANCE  RATI 0( MU)  -  0.215 

INDUCED  VELCC I  TV  ( siu )  -  16.51  -~/S 

INFLOW  RAT  IQ (LAMBDA)  -  -0.151c 

FLAT  PLATE  DRAG  -  272 . 0  LB 

THRUST  -  9607.9  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.01016 

BLADE  SOLIDITY  -  0.0708 

COLLECTIVE  (THETA  0)  -  17.9  DEG 

CONING  ANGLE  (aO)  -  13.6  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (el)  -  6.0  DEG 


0.01016 


6.0  DEG 


LATERAL  FLAPPING  ANGLE  (bl)  -  3.9  DEG 

DISC  ANGLE  OF  ATTACK  -  -7.6  DEG 

INDUCED  POWER  -  299.3  HP 

PARASITE  POWER  -  76 . 8  HP 

PROFILE  POWER  -  195.8  HP 

TOTAL  POWER  -  560.9  HP 

LOCAL  ANGLE  OF  ATTACK ( 027ODEG&4OS  RADIUS  )  -  11.6  DEG 

LOCAL  ANGLE  OF  ATTACK(027ODEG&TIP)  -  13.7  DEG 

LOCAL  ANGLE  OF  ATTACK  (027ODEG&PT . OF  UT-0.4)  -  14.7  DEG 

RADIAL  POINT  OF  UT-0.4  (@270DEG)  -  0.615 

1  RESULTS  FOR  RANGE  OF  334.8  MILES,  ARE  AS  FOLLOWS  : 

AIRSPEED  -  142.00  KNOTS 

ADVANCE  RATIO(MU)  -  0.325 

INDUCED  VELOCITY(NU)  -  10.80  FT/S 

INFLOW  RAT I 0( LAMBDA)  -  -0.0930 

FLAT  PLATE  DRAG  -  647.9  LB 

THRUST  -  9625.0  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.01018 

BLADE  SOLIDITY  -  0.0788 

COLLECTIVE  (THETA  0)  -  19.9  DEG 

CONING  ANGLE  (aO)  -  13.4  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (al)  -  9.6  DEG 

LATERAL  FLAPPING  ANGLE  (bl)  -  5.7  DEG 

DISC  ANGLE  OF  ATTACK  -  -13.5  DEG 

INDUCED  POWER  -  189.0  HP 

PARASITE  POWER  -  282.3  HP 

PROFILE  POWER  -  215.7  HP 

TOTAL  POWER  -  687.0  HP 

LOCAL  ANGLE  OF  ATTACK (027ODEG&4OS  RAOIUS  )  -  -2.2  DEG 

LOCAL  ANGLE  OF  ATTACK (027ODEG&TIP)  -  18.3  DEG 

LOCAL  ANGLE  OF  ATTACK  (027ODEG&PT. OF  UT-0.4)  -  18.2  DEG 
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From  these  results,  it  can  be  seen  that  the  velocity  for 
m-i niimm  power  is  92  knots  (155.3  ft/s)  . 

From  equation  (23)  ,  y  ■  arcsin  (V  /V) 

C  V 

»  arcsin  (1.67/155.3) 

-  0.615° 

'POLAR*  is  then  run  as  shown  below: 


.RU  POLAR 


TITLE  (TWO  LINES  OF  UP  TO  60  CHARACTERS) 

:  SERUIUE  CEILING  ESTIMATE 

ARE  UNITS  IN  IMPERIAL  OR  METRIC  (I  OR  M)  ?  I 

ALL  UP  WEIGHT  (N  OR  LS)  -  ?  9604 

ROTOR  TIP  SPEED  (M/S  OR  FT/S)  -  ?  717.2 

ROTOR  RAO  I US  (M  OR  FT)  -  7  21 

ROTOR  BLADE  CHORD  (M  OR  FT)  -  ?  1.3 

NUMBER  OF  ROTOR  BLADES  -  ?  4 

2D  LIFT  CURUE  SLOPE  (/RAO)  -  ?  5.73 

TIP  LOSS  FACTOR  -  ?  0.97 

SLADE  TWIST  (DEG)  -  ?  -8 

EQUIUALENT  FLAT  PLATE  AREA  <M**2  OR  FT**2 )  -  ?  17 
LOCKS  NUMBER  -  ?  15 

DRAG  POLAR  CO-EFFICIENT  (DELO)  -  ?  0.0087 
DELTA1  (/RAD)  -  ?  -0.0216 
DELTA2  (/RAD**2)  -  ?  0.4 

AIR  DENSITY  (KG/M**3  OR  SLUG/FT**3 )  -  ?  0.001327 

IS  RANGE/ENOURANCE  REQUIRED  ?  (Y  OR  N)  :  N 

IS  A  SPEED-POWER  POLAR  REQUIRED  ?  (Y  OR  N)  :  N 

IS  POWER  REQUIRED  FOR  A  GIUEN  CLIMB  ANGLE 
i.e.  RATE  OF  CLIMB  (Y  OR  N)  ?  Y 

VELOCITY  IN  KNOTS  ?  s  92 

WHAT  IS  CLIMB  ANGLE  -  GAMMA  -  ASIN(VU/U)  -  (DEG)  ?  :  0.615 
STOP 

END  OF  EXECUTION 

CPU  TIME:  1.38  ELAPSED  TIME:  53.96 
EXIT 
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.TV  POLAR. OUT 

§ERUIUE  CEILING  ESTIMATE 

ALL  UP  WEIGHT  -  9604.0  LB 

AIR  DENSITY  -  0.001327  SLUG/FT**3 

TIP  SPEED  -  717.2  FT/S 

ROTOR  RADIUS  -  21.0  FT 

ROTOR  CHORD  -  1 . 3  FT 

NUMBER  OF  ROTOR  BLADES  -  4.0 

LIFT  CURUE  SLOPE  -  5.73  /RAO 

TIP  LOSS  FACTOR  -  0.97 

BLADE  TWIST  -  -8.0  DEG 

EQUIUALENT  FLAT  PLATE  AREA  -  17.0  FT**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

DELTA  1  -  -0.0216  /RAD 
DELTA  2  -  0.4000  /RAO**2 

AIRSPEED  -  92.00  KNOTS 

ADUANCE  RAT  IOC  MU )  -  0.214 

INDUCED  UELOCITYCNU)  -  16.49  FT/S 

INFLOW  RAT  IOC  LAMBDA )  -  -0.0541 

FLAT  PLATE  DRAG  -  272.0  LB 

THRUST  -  9610.8  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.01016 

BLADE  SOLIDITY  -  0.0788. 

COLLECT  I UE  (THETA  0)  -  18.1  DEG 

CONING  ANGLE  CaO)  -  13.6  DEG 

LONGITUDINAL  FLAPPING  ANGLE  Cal)  -  6.0  DEG 

LATERAL  FLAPPING  ANGLE  (bl)  -  3.9  DEG 

DISC  ANGLE  OF  ATTACK  -  -8.2  DEG 

INDUCED  POWER  -  298.2  HP 

parasite  power  -  76.3  hp 

PROFILE  =CWER  -  196.  2  -P 

TOTAL  POWER  -  561.1  HP 

REQUIRED/ AVAILABLE  SHAFT  POWER  -  590.2  HP 

CLIMB  POWER  -  29.1  HP 

RATE  OF  CLIMB  -  100.0  FT/MIN 

CLIMB  ANGLE  -  0.6  DEG 

LOCAL  ANGLE  OF  ATTACK ( Q2700EG&40J<  RADIUS  )  -  1 

LOCAL  ANGLE  OF  ATTACK  C  [3270DEG&T  IP  )  -  13.8  DEG 

LOCAL  ANGLE  OF  ATTACK  X 0270DEGSPT . OF  UT-0.4)  - 
RAO  I AL  POINT  OF  UT-0.4  (@270DEG)  -  0.614 


The  relevant  parameters  for  this  case  are: 


P,  -  288.2  Hp,  P  -  76.8  Hp,  P  -  29.1  Hp, 
l.  p  c 


?°calc  "  196’2  HP>  °<1»270>  "  13'8° 


•*.P  -  196.2  (1  +  13.-..a-.-~  -12  )  .  284.5  Hp 

o  4 


1.2  DEG 
14.7  DEG 


RU  POLAR 


TITLE  (TWO  LINES  OF  UP  TO  60  CHARACTERS) 

:  SLADE  ELEMENT  METHOD  -  IMPERIAL  UNITS 
HOUER  CEILING  O.G.E. 

ARE  UNITS  IN  IMPERIAL  OR  METRIC  ( I  OR  M)  ?  I 

ALL  UP  WEIGHT  (N  OR  L8)  -  ?  9604 

ROTOR  TIP  SPEED  (M/S  OR  FT/S)  -  ?  717.2 

ROTOR  RADIUS  CM  OR  FT)  -  ?  21 

ROTOR  BLADE  CHORD  (M  OR  FT)  -  ?  1.3 

NUMBER  OF  ROTOR  SLADES  -  7  4 

2D  LIFT  CURUE  SLOPE  (/RAD)  -  7  5.73 

TIP  LOSS  FACTOR  -  7  0,97 

SLAOE  TWIST  (DEG)  -  ?  -3 

ECU  I UALENT  FLAT  PLATE  AREA  ( OR  F7*C  >  -  7  17 
LOCKS  NUMBER  -  ?  15 

DRAG  POLAR  CO-EFFICIENT  (DELO)  -  ?  0.0087 
DELTA 1  (/RAO)  -  ?  -0.0216 
DELTA2  (/RAD**2)  -  ?  0.4 

AIR  DENSITY  (KG/M**3  OR  SLUG/FT**3 )  -  ?  0.002041 

IS  RANGE/ENDURANCE  REQUIRED  ?  (Y  OR  N)  :  N 

IS  A  SPEED-POWER  POLAR  REQUIRED  ?  (Y  OR  N)  :  N 

IS  POWER  REQUIRED  FOR  A  GIUEN  CLIMB  ANGLE 
i.e.  RATE  OF  CLIMB  (Y  OR  N)  ?  Y 

UELOCITY  IN  KNOTS  ?  :  0 

WHAT  IS  CLIMB  ANGLE  -  GAMMA  -  ASINCW/U)  -  (DEG)  ? 
STOP 

END  OF  EXECUTION 

CPU  TIME:  1.15  ELAPSED  TIME:  37.68 
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BLADE  ELEMENT  METHOD  -  IMPERIAL  UNITS 
HOUER  CEILING  O.G.E. 

ALL  UP  WEIGHT  -  9604.0  LB 

AIR  DENSITY  -  0.002041  SLUG/FT**3 

TIP  SPEED  -  717.2  FT/'S 

ROTOR  RADIUS  -  21.0  FT 

ROTOR  CHORD  -  1.3  FT 

NUMBER  OF  ROTOR  BLADES  -  4.0 

LIFT  CURUE  SLOPE  -  5.73  /RAD 

TIP  LOSS  FACTOR  -  0.97 

BLADE  TWIST  -  -0.0  DEG 

EQUIUALENT  FLAT  PLATE  AREA  -  17.0  FT**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

DELTA  1  -  -0.0216  /RAO 
DELTA  2  -  0.4000  /RAD**2 

AIRSPEED  -  0.00  KNOTS 

ADUANCE  RATI 0( MU)  -  0.000 

INDUCED  UELOCITY(NU)  -  53.56  FT/S 

INFLOW  RATI 0( LAMBDA)  -  -0.0747 

FLAT  PLATE  DRAG  -  0 . 0  LB 

THRUST  -  9604.0  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.00660 

BLADE  SOLIDITY  -  0.0780 

COLLECT  I UE  (THETA  0)  -  17.9  DEG 

CONING  ANGLE  (aO)  -  9.6  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (al)  -  0.0  DEG 

LATERAL  FLAPPING  ANGLE  (bl)  -  0.0  DEG 

DISC  ANGLE  OF  ATTACK  -  0.0  DEG 

INDUCED  POWER  -  935.2  HP 

3ARASI”E  rower  «  0.0  UP 

PROFILE  POWER  »  204.3  HP 

total  POWER  -  1139.4  HP 

REQU I RED/AU A I LASLE  SHAFT  POWER  -  1139.4  HP 

CLIMB  POWER  -  0.0L  HP 

RATE  OF  CLIMB  -  0.0  FT/MIN 

CLIMB  ANGLE  -  0.0  DEG 

LOCAL  ANGLE  OF  ATTACK ( 027ODEG&4QS  RADIUS  )  - 
LOCAL  ANGLE  OF  ATTACK ( @270 DEG&T IP )  -  5.6  DEG 

LOCAL  ANGLE  OF  ATTACK  ( 027ODEG&PT. OF  UT-0.4)  - 
RADIAL  POINT  OF  UT-0.4  (@270DEG)  -  0.400 


4.1  DEG 
4.1  DEG 


Available  shaft  power  is  given  by: 

Pavail  *  1275  (1'13  *  0*002041/0.002259  -  0.13) 

■  1139  Hp. 

Thus  hover  ceiling  O.G.E.  is  approximately  3400  ft. 


4.4.4.  Hover  ceiling  In  Ground  Effect 


-  I 

A  ground  effect  relationship  (ref.  8)  is  given  as: 


1 


/R_xa 

ja|lj  MZ' 

4Ct  l+(y/A)a 

For  example,  taking  z  =  R/2  and  using  the  hover  O.G.E.  case 

as  a  basis: 

T_  _  _ 1 _ 

*  0.0788  x  5.73  x  0.0747  (1/2) 2 

1  4  x  0.00660  1 


*  1.4693. 


A  helicopter  of  weight  9604  lbs  hovering  I.G.E.  is  therefore 
equivalent  to  a  helicopter  of  weight  9604/1.4693  *  6536  lbs  hovering 
O.G.E. 


'POLAR'  is  run  as  for  §4.4.3  at  this  equivalent  weight  for 
an  altitude  of  14000  ft,  and  is  the  final  result  after  several  iterations 


.TV  POLAR. OUT 

Slade  element  method  -  imperial  units 

HOUER  CEILING  I.G.E. 

ALL  UP  WEIGHT  -  6536.0  LB 

AIR  DENSITY  -  0.001461  SLUG/FT**3 

TIP  SPEED  -  717.2  FT/S 

ROTOR  RAO  I US  -  21.0  FT 

ROTOR  CHORD  -  1 . 3  FT 

NUMBER  OF  ROTOR  BLADES  -  4.0 

LIFT  CURUE  SLOPE  -  5.73  /RAO 

TIP  LOSS  FACTOR  -  0.97 

BLADE  TWIST  -  -8.0  DEG 

EQU I UALENT  FLAT  PLATE  AREA  -  17.0  FT**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

DELTA  1  -  -0.0216  /RAD 
DELTA  2  -  0.4000  /RAO**2 

AIRSPEED  -  0.00  KNOTS 

ADUANCE  RATIO! MU)  -  0.000 

INOUCED  UELOCITY(NU)  -  52.22  FT/S 

INFLOW  RAT IO (LAMBDA)  -  -0.0728 

FLAT  PLATE  DRAG  -  0 . 0  LB 

THRUST  -  6536.0  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.00628 

BLADE  SOLIDITY  -  0.0788 

COLLECTIVE  (THETA  0)  -  17.5  DEG 

CONING  ANGLE  (aO)  -  9.1  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (al)  -  0.0  DEG 

LATERAL  FLAPPING  ANGLE  (bl)  -  0.0  DEG 

DISC  ANGLE  OF  ATTACK  -  0.0  DEG 

INDUCED  POWER  -  620:6  HP 

PARASITE  POWER  -  0 . 0  HP 

PROFILE  POWER  -  142.1  HF 

TOTAL  POWER  -  762.7  HP 

REQUIRED/ AVAILABLE  SHAF-  POWER  -  762.7  HP 

CLIMB  POWER  -  0.0  HP 

RATE  OF  CLIMB  -  0.0  FT/M  IN 

CLIMB  ANGLE  -  0.0  DEG 

LOCAL  ANGLE  OF  ATTACK ( (a270DEG&4GS  RADIUS  )  - 
LOCAL  ANGLE  OF  ATTACK ( 027ODEG&T I P )  -  5.3  DEG 

LOCAL  ANGLE  OF  ATTACK  (027ODEG&PT . OF  UT-0.4)  - 
RADIAL  POINT  OF  UT-0.4  I@270OEG)  -  0.400 


3.9  DEG 
3.9  DEG 


Here  P  »  763  Hp. 
req 

At  14,000  ft, 

Pavail  "  1275  {1*13  x  0-001461/0.002259  -  0.13) 

-  766  Hp. 

Thus  the  Hover  Ceiling  I.G.E.  is  approximately  14,000  ft. 

It  should  be  noted  that  in  some  situations,  a  helicopter  can 
hover  I.G.E.,  but  has  insufficient  power  to  move  forward  without 
decending.  A  check  to  see  if  this  occurs  can  be  made  in  the  following 
manner: 


(1) 


Take  a  case  for  low  u  (say  0.02)  .in  level  flight  at  the 
hover  O.G.E.  ceiling. 


(2) 


Calculate  the  I.G.E.  thrust  ratio 


(3)  Calculate  the  equivalent  O.G.E.  weight. 

(4)  Run  'POLAR'  for  that  p  at  the  hover  I.G.E.  ceiling. 

(5)  If  the  power  required  is  greater  than  the  power  available 
then  the  helicopter  cam  hover  but  cannot  move  forward 
without  descending. 

(6)  Repeat  steps  1  5  over  a  range  of  y,  if  required. 

An  example  showing  the  results  of  this  procedure  is  set  out 

below: 


SLADE^LEM^NT  METHOD  -  IMPERIAL  UNITS 
MU  -  0.02,  O.G.E. 

ALL  UP  WEIGHT  -  9604.0  LB 

AIR  DENSITY  -  0.002041  SLUG/FT**3 

TIP  SPEED  -  717.2  FT/S 

ROTOR  RAO  I US  -  21.0  FT 

ROTOR  CHORD  -  1.3  FT 

DUMBER  OF  ROTOR  SLADES  -  4.0 

;_IFT  CURLE  SLOPE  *  r . 73  /FAQ 

TIP  LOSS  -ACTOR  -  0.97 

BLADE  TWIST  -  -8.0  DEG 

EGUIUALENT  FLAT  PLATE  AREA  -  17.0  FT**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

DELTA  1  -  -0.0216  /RAO 
DELTA  2  -  0.4000  /RAO**2 

AIRSPEED  -  8.50  KNOTS 


17  FT/S 

0701 

LB 

0 . 00660 


0.5 

DEG 


ADUANCE  RAT  IQ ( MU )  -  0.020 

INDUCED  UELGCITY(NU)  -  50.17  FT/S 

INFLOW  RAT  IOC  LAMBDA )  -  -0.0701 

FLAT  PLATE  DRAG  -  3.6  LB 

THRUST  -  9604.0  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.00660 

BLADE  SOLIDITY  -  0.0738 

COLLECT  I UE  (THETA  0)  -  17.5  DEG 

CONING  ANGLE  <a0)  -  9.6  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (al)  -  0.5 

LATERAL  FLAPPING  ANGLE  (bl)  -  0.3  DEG 

DISC  ANGLE  OF  ATTACK  -  -0.5  DEG 

INDUCED  POWER  -  876.1  HP 

PARASITE  POWER  -  0 . 1  HP 

PROFILE  POWER  -  202.9  HP 

TOTAL  POWER  -  1079.1  HP 

REQU I RED/AUA I LABLE  SHAFT  POWER  -  1079.1 

CLIMB  POWER  -  0.0  HP 

RATE  OF  CL  I M8  -  0.0  FT/M  IN 

CLIMB  ANGLE  -  0.0  DEG 

LOCAL -ANGLE  OF  ATTACK 082 7 0DEG&40H  RADIUS 

LOCAL  ANGLE  OF  ATTACK ( Q270DEG&T I P )  -  5 

LOCAL  ANGLE  OF  ATTACK  ( 027ODEG&PT . OF  UT- 

RAOIAL  POINT  OF  UT-0.4  T027ODEG)  -  0.42 


270DEG&T IP)  -  5.9 
027ODEG&PT . OF  UT-0. 
(027ODEG)  -  0.420 


S  )  - 
5.9  DEG 
*0.4)  - 


Calculation  of  ground  effect  factor  gives: 


1  - 


0.0788  x  5.73  x  0.0701 
4  x  0.0066 


(1/2) 


1.3835 


l+(0. 02/0. 0701) 


Thus  equivalent  O.G.E.  weight  is  9604/1.3835  *  6942  lbs. 


BLADE  ELEMENT  METHOD  -  IMPERIAL  UNITS 
MU  -  0.02.  I.G.E. 

ALL  UP  WEIGHT  -  6942.0  LB 

AIR  DENSITY  -  0.001461  SLUG/FT<**3 

TIP  SPEED  -  717.2  FT/S 

ROTOR  RADIUS  -  21.0  FT 

ROTOR  CHORD  -  1 . 3  FT 

NUMBER  OF  ROTOR  BLADES  -  4.0 

LIFT  CURUE  SLOPE  -  5.73  /RAD 

TIP  LOSS  FACTOR  -  0.97 

BLADE  TWIST  -  -8.0  DEG 

EQUIUALENT  FLAT  PLATE  AREA  -  17.0  FT**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

DELTA  1  -  -0.0216  /RAO 
DELTA  2  -  0.4000  /RAO**2 

AIRSPEED  -  8.50  KNOTS 

ADUANCE  PAT  IQ (MU)  -  0.020 

I NOUCSD  oELOC I  TV ( NU )  -  5  0.^5  FT/S 

INFLOW  RAT  1 0( LAMBDA)  -  -0.0705 

FLAT  PLATE  ORAG  -  2 . 6  LS 

THRUST  -  6942.0  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.00667 

BLADE  SOLIDITY  -  0.0788 

COLLECT  I UE  (THETA  0)  -  17.6  DEG 

CONING  ANGLE  (aO)  -  9.6  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (al)  -  0.5  DEG 

LATERAL  FLAPPING  ANGLE  (bl)  -  0.3  DEG 

DISC  ANGLE  OF  ATTACK  -  -0.5  DEG 

INDUCED  POWER  -  636.6  HP 

PARASITE  POWER  -  0 . 1  HP 

PROFILE  POWER  -  146.0  HP 

TOTAL  POWER  -  782.7  HP 

REQU I RED/AUA I LABLE  SHAFT  POWER  -  782.7  HP 

CLIMB  POWER  -  0.0  HP 

RATE  OF  CLIMB  -  0.0  FT/MIN 

CLIMB  ANGLE  -  0.0  DEG 

LOCAL  ANGLE  OF  ATTACK ( 027OOEG&4OH  RADIUS  )  - 
LOCAL  ANGLE  OF  ATTACK(§270OEG&TIP)  -  5.9  DEG 

LOCAL  ANGLE  OF  ATTACK  ( 027ODEG&PT . OF  UT-0.4)  - 
RAO  I AL  POINT  OF  UT-0.4  C@270DEG)  -  0.420 


4.3  DEG 
4.7  DEG 


Now  P  is  783  Hp,  but  available  power  at  14,000  ft  is 
req 

766  Hp  (found  previously) . 


Thus  if  helicopter  moves  forward,  it  will  descend. 
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5.  COMPARISON  OF  'POLAR'  AND  OTHER  METHODS  OF  PERFORMANCE  ESTIMATION 


As  mentioned  in  Section  1,  the  methods  of  performance 
estimation  available  prior  to  'POLAR'  were: 

1.  Program  'ENERGY'  (see  Appendix  4),  derived  from 
'MCEP',  and 

2.  Performance  Charts  of  References  1  and  2. 

This  section  compares  results  obtained  for  the  hypothetical  helicopter 
of  Section  4,  using  'POLAR'  and  the  above-mentioned  methods. 

5.1  'POLAR'  and  'ENERGY' 

Figure  4  shows  an  overlay  plot  of  speed  against  power,  with 
a  breakdown  of  power  components,  using  'POLAR'  and  'ENERGY'. 

As  can  be  seen  from  this  figure,  agreement  of  the  power 
components  is  generally  good.  The  main  discrepancy  between  'ENERGY' 
and  'POLAR'  occurs  at  high  speed.  This  is  partially  due  to  'POLAR' 
predicting  a  lower  profile  power,  but  mostly  due  to  compressibility 
power  calculated  by  'ENERGY'.  'POLAR'  does  not  take  compressibility 
effects  into  account. 

5.2  'POLAR*  and  Performance  Charts 

Figures  5,  6  and  7  show  power-speed  polars  obtained  using 
'POLAR',  together  with  points  obtained  from  the  performance  charts  of 
references  1  and  2,  for  rotor  blade  twists  of  0,  8  and  16  degrees 
respectively . 

5.2.1  'POLAR'  and  Performance  Charts  of  Gessow  and  Tapscott 

From  Figure  5,  it  can  be  seen  that  agreement  is  good  between 
'POLAR'  and  the  Reference  2  charts,  except  at  low  speed.  This  is  due  to 
the  fact  that  'POLAR'  applies  a  semi -empirical  factor  to  the  induced 
velocity  (for  uS0.14)  whereas  the  charts  use  only  the  momentum  theory 
value  (see  §2.2).  Thus  the  charts  always  underestimate  the  induced 
velocity  (induced  power)  at  low  speed. 


Figures  6  and  7  also  show  good  agreement,  except  at  low  speed 
where  the  induced  power  is  underestimated  by  the  charts,  as  expected. 

For  both  cases  (8  and  16°  twist),  at  high  speed  'POLAR*  underestimates 
the  profile  power  because  'POLAR'  ignores  the  high  order  terms  in 
U  and  8  ,  in  the  profile  power  expression  (see  §2.5).  Thus  at  high 
speeds  (y)»  neglect  of  these  terms  comes  into  effect.  For  example,  at 
Vi  -  0.3,  these  terms  cause  a  13%  difference  in  profile  power  for  8  and 
16  degrees  twist,  but  this  represents  only  a  5%  difference  in  total  power 

5.2.2  'POLAR'  and  Performance  Charts  of  Tanner 

The  performance  charts  of  Reference  1  are  applicable  only  to 
hovering  helicopters  and  aircraft  flying  at  p£0.25  for  8°  rotor  twist 
and  pS0.3  for  0°  twist. 

Figure  5  shows  poor  agreement  for  the  total  power.  In  the  hover 
case  this  is  due  to  the  semi-empirical  factor  applied  at  low  speed  by 
'POLAR' (as  in  §5.2.1).  At  p  =  0.3,  the  charts  of  Reference  1  indicate 
that  the  rotor  is  stalled.  This  discrepancy  is  thus  probably  due  to 
Tanner's  inclusion  of  compressibility  and  stall  effects,  which  are 
neglected  in  the  case  of  'POLAR'. 

Figure  6  again  shows  poor  agreement  at  hover,  for  the  reason 
mentioned  previously.  However,  at  p  =  0.25  and  p  =  0.3  we  can  see 
excellent  agreement.  In  this  case  the  rotor  is  not  stalled,  thus 
verifying  that  the  difference  in  Figure  5  is  due  to  Tanner's  stall 
effect. 

6.  CONCLUDING  REMARKS 

Program  'POLAR'  has  been  developed  to  give  a  wide  range  of 
performance  calculations  within  the  operational  boundaries  of  most 
conventional  helicopters. 

'POLAR'  allows  rapid  estimation  and  is  more  generally 
applicable  than  existing  performance  charts. 

Being  based  on  blade  element  theory,  'POLAR'  overcomes  the 
deficiencies  of  the  energy  balance  method  (program  'ENERGY')  in  that 
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blade  twist  is  taken  into  account  and  blade  stall  is  indicated  by 
calculation  of  blade  angles  of  attack  at  critical  locations  on  the 
rotor  disk. 

Since  'POLAR'  does  not  take  compressibility  and  stall  effects 
into  account,  it  is  limited  to  y£0.3. 
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APPENDIX  1  -  USING  PROGRAM  'POLAR' 


I 


i 


p 


Al.l  Loading  'POLAR'  on  the  DEC-10 

For  ARL  users,  a  description  is  shown  below  of  how  to 
load  program  ’POLAR'  on  the  DEC-10  computer.  'POLAR'  is  also  available 
on  the  VAX-11  computer  and  is  loaded  in  a  similar  fashion. 


.PLE  MX  D/039 

MESSAGE  TRANSMITTED  at  13:29 


j:OPR:  -  D/039  MX  DTAO 
.AS  DTAO 


DTAOOO  assigned 


.DIR  DTAO: 

Tape  ID:  D/039 
Free:  423  biks, 
POLAR  . CCL  1 
POLAR  .FOR  60 
DING  .REL  1 
ENERGY . CCL  1 
ENERGY. FOR  32 
SLADE  . “CP  -3 
TSUea 1 . REL  1 0 
CPU  . REL  1 
ATMOS  .FOR  15 


13  files 
2Q-Apr-83 
7-Sep-83 
14- Apr-33 
20-Apr-33 
7-Seo-83 
a-:un-33 
3  0-Mar-33 
9-Nov-79 
23- Jun-83 


. COPY/X-OTAO : POLAR . * , TSUB8 1 . * , CPU . * , D I NG . * 

. COM  POLAR 

FORTRAN:  POLAR 

POLAR 

ELMENT 

LABDA 

TAN 

IMPLIC 

COEFF 

POWER 

ALPH 

ROC 

CLIMB 

RNGEND 

.TY  POLAR. CCL 

POLAR , TSUB8 1 , CPU .DING 

POLAR/SAU- 

/GO 

.R  LINK 
♦@POLAR 
EXIT 

.DIR  POLAR . EXE 
POLAR  EXE  60  <057> 


7-SEP-83 


DSKD:  11031,1063] 


When  running  'POLAR1 ,  the  results  are  output  onto  one  of 
two  files.  If  performance  is  required  at  one  particular  airspeed, 
or  range/endurance  cure  desired,  the  results  are  output  onto  file 
'POLAR. OUT'.  If  a  speed-power  polar  is  required,  the  results  are  output 
onto  file  'POLAR.DAT' ,  which  is  then  input  to  program  'TRANS'  for 
presentation  in  column  or  graphical  form  with  velocity  as  the  independent 
variable.  An  example  of  this  procedure  is  given  in  Appendix  2. 

A1.3  Renaming  Files 

A  file  renaming  feature  is  included  in  'POLAR1  to  allow 
up  to  five  consecutive  runs  to  be  made  without  overwriting  the 
results  obtained  from  a  previous  run.  As  stated  above,  results  for 
each  run  are  always  output  to  files  POLAR. OUT  or  POLAR.DAT  If  a  second 
run  is  made,  the  preceeding  POLAR. OUT  or  POLAR.DAT  will  automatically 
be  renamed  P0L1.0UT  or  POLl.DAT.  This  process  will  continue  up  to 
P0L4.0UT  or  POL4.DAT. 


A1.4  Special  Cases 

It  was  found  that  for  some  flight  conditions,  when 
calculating  rate  of  climb  for  a  given  shaft  power  and  airspeed, 
convergence  would  not  occur  unless  the  climb  angle,  could  be  roughly 
estimated  initially.  Therefore,  'POLAR'  calculates  the  power  required 
in  level  flight  411(1  the  powex  required  in  a  vertical  climb 

(P  and  uses  a  linear  interpolation  to  initially  approximate  yc» 


i.e. 


Y  =  (P  ..-P,  ,)/(P  -P  90° 

c  avail  level  vert  avail 


In  some  cases  when  the  rate  of  climb  is  required  for  a 
given  airspeed  and  shaft  power,  the  power  may  not  be  sufficient  to 
sustain  a  climb  and  therefore  the  helicopter  descends.  'POLAR' 
outputs  a  message  to  this  effect  which  is  demonstrated  below. 


.RU  POLAR 


TITLE  (TWO  LINES  OF  UP  TO  60  CHARACTERS) 
:  SLADE  ELEMENT  METHOD  -  IMPERIAL  UNITS 
DESCENDING  FLIGHT 


?  I 


ARE  UNITS  IN  IMPERIAL  OR  METRIC  (I  OR  M) 

ALL  UP  WEIGHT  (N  OR  LB)  -  ?  9604 
ROTOR  TIP  SPEED  (M/S  OR  FT/S)  -  ?  717.2 
ROTOR  RADIUS  (M  OR  FT)  -  ?  21 
ROTOR  SLADE  CHORD  (M  OR  FT)  -  ?  1.3 
NUMBER  OF  ROTOR  BLADES  -  ?  4 
2D  LIFT  CURUE  SLOPE  (/RAO)  -  ?  5.73 
TIP  LOSS  FACTOR  -  ?  0.97 
BLADE  TWIST  (DEG)  «  ?  -3 

EQUIVALENT  FLAT  PLATE  AREA  ■' M**2  OR  FT**2 ) 

LOCKS  NUMBER  -  ?  1? 

DRAG  POLAR  CO-EFFICIENT  (DELO)  -  ?  0.0087 
DELTA 1  (/RAO)  -  ?  -0.0216 
DELTA2  (/RAD*+2)  -  ?  0.4 

AIR  DENSITY  (KG/M**3  OR  SLUG/FT**3 )  -  ?  0.002259 
IS  RANGE/ENOURANCE  REQUIRED  ?  (Y  OR  N)  :  N 
IS  A  SPEED- POWER  POLAR  REQUIRED  ?  (Y  OR  N)  :  N 


?  17 


IS  POWER  REQUIRED  FOR  A  GIVEN  CLIMB  ANGLE 
i.e.  RATE  OF  CLIMB  (Y  OR  N)  ?  N 


VELOCITY  IN  KNOTS  ?  :  12 

WHAT  IS  POWER  AVAILABLE  AT  MAIN  ROTOR  SHAFT  (HP)  ? 


STOP 


ENO  OF  EXECUTION 

CPU  TIME !  1.27  ELAPSED  TIME:  36.78 


' HEL?8hpT&SUHAS  INSUFFICIENT  POWER  FOR  LEUEL  FLIGHT 
i.a.  AIRCRAFT  IS  DESCENDING 

BLADE  ELEMENT  METHOD  -  IMPERIAL  UNITS 
DESCENDING  FLIGHT 
ALL  UP  LIE  I GHT  -  96  04.0  LB 

AIR  DENSITY  -  0.002259  SLUG/FT**3 

TIP  SPEED  -  717.2  FT/'S 

ROTOR  RADIUS  -  21.0  FT 

ROTOR  CHORD  -  1.3  FT 

NUMBER  OF  ROTOR  BLADES  «  4.0 

LIFT  CURUE  SLOPE  -  5.73  /RAD 

TIP  LOSS  FACTOR  -  0.97 

BLADE  TWIST  -  -8.0  DEG 

EQUIUALENT  FLAT  PLATE  AREA  -  17.0  FJ**2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

DELTA  1  -  -0.0216  /RAD 
DELTA  2  -  0.4000  /RAD**2 

AIRSPEED  -  12.00  KNOTS 

ADUANCE  RAT IO (MU)  -  0.028 

INDUCED  UELOCITY(NU)  -  45.48  FT/S 

INFLOW  RAT  I O ( LAMBDA )  -  -0.0632 

FLAT  PLATE  DRAG  -  7.9  LB 

THRUST  -  9603.9  LB 

THRUST  CO-EFFICIENT  (CT)  -  0.00597 

BLAOE  SOLIDITY  -  0.0788 

COLLECT  I UE  (THETA  0)  -  16.4  DEG 

CONING  ANGLE  (aO)  -  9.6  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (al)  -  0.6  DEG 


LATERAL  FLAPPING  ANGLE  (bl)  -  0.3  DEG 

D ; SC  ANGLE  OF  A i . ACX  —  0.4  CEG 

INDUCED  POWER  -  794.2  HP 

PARASITE  POWER  -  0.3  HP 

PROFILE  POWER  -  211.9  HP 

TOTAL  POWER  -  1006.4  HP 

REQU I RED/AUA I LA8LE  SHAFT  POWER  -  1000.0  HP 

CLIMB  POWER  «  -6.4  HP 

RATE  OF  CLIMB  -  -22.1  FT/M IN 

CLIM8  ANGLE  -  -1.0  DEG 

LOCAL  ANGLE  OF  ATTACK ( 027OOEG&4OX  RADIUS  )  -  4.1  DEG 

LOCAL  ANGLE  OF  ATTACK ( 027OOEG4T IP )  -  5.2  DEG 

LOCAL  ANGLE  OF  ATTACK  ( 027ODEG&PT . OF  UT-0.4)  -  4.6  DEG 

RADIAL  POINT  OF  UT-0.4  C027ODEG)  -  0.428 


A  further  case  to  be  considered  is  when  the  rate  of  climb  is 
required  for  a  given  airspeed  and  shaft  power,  and  the  available  power  is 
greater  than  that  required  for  vertical  flight.  'POLAR'  then  outputs  a 
message  stating  the  excess  power  available  and  then  gives  results  for  a 
vertical  climb  at  the  specified  airspeed.  In  the  example  below, 
the  inputs  are  identical  to  the  previous  example  except  that  the  power 
available  is  now  1275  Hp.  Only  the  output  file  is  given. 


'HELI  2hpTE8UHAS  12.9  HP  MORE  THAN  IS 
REQUIRED  FOR  A  UERTICAL  CL  I MS.  RESULTS  BELOW 
ARE  THOSE  FOR  A  UERTICAL  CLIMB. 

BLADE  ELEMENT  METHOO  -  IMPERIAL  UNITS 
EXCESS  POWER  DEMONSTRATION 
ALL  UP  WEIGHT  -  9604.0  LB 

AIR  DENSITY  -  0.002259  SLUG/FT**3 

TIP  SPEED  -  717.2  FT/S 

ROTOR  RADIUS  -  21.0  FT 

ROTOR  CHORD  -  1 . 3  FT 

NUMBER  OF  ROTOR  BLADES  -  4.0 

LIFT  CURUE  SLOPE  -  5.73  /RAO 

TIP  LOSS  FACTOR  -  0.97 

BLADE  TWIST  -  -8.0  DEG 

ECU IUALENT  FLAT  PLATE  AREA  -  17.0  "«2 

LOCKS  NUMBER  -  15.0 

DRAG  POLAR  CO-EFFICIENT  (DELTA  0)  -  0.0087 

OELTA  1  -  -0.0216  /RAO 
DELTA  2  -  0.4000  /RAO**2 

AIRSPEED  -  12.00  KNOTS 


AOUANCE  RAT  IOC  MU)  - 
INDUCED  UELOCITY(NU) 

INFLOW  RATIOCLAMBDA) 

FLAT  PLATE  DRAG  - 
THRUST  -  9612.4  LB 

THRUST  CO-EFFICIENT  (CT)  - 
BLADE  SOLIDITY  -  0.0788 

COLLECT  I UE  (THETA  0)  -  18 

CONING  ANGLE  (aO)  -  8.9  DEG 

LONGITUDINAL  FLAPPING  ANGLE  (el) 


000 

39 . 44  FT/S 
-0.0832 
7.9  LB 

0.00597 


1  DEG 


0.0  DEG 


LATERAL  FLAPPING  ANGLE  <bl)  -  0.0  DEG 

DISC  ANGLE  OF  ATTACK  -  -90.0  DEG 
INDUCED  POWER  -  689.3  HP 

PARASITE  POWER  -  0 . 3  HP 

PROFILE  POWER  -  218.9  HP 

TOTAL  POWER  -  908.5  HP 

REQU I RED/AUA I LABLE  SHAFT  POWER  -  1275.0  HP 

CLIMB  POWER  -  353.7  HP 

RATE  OF  CLIM8  -  1215.2  FT/MIN 

CLIMB  ANGLE  -  90.0  DEG 

LOCAL  ANGLE  OF  ATTACK  ( (32  70DEG&4  OH  RADIUS  )  - 
LOCAL  ANGLE  OF  ATTACK(fi270DEG&TIP)  ■  5.3  DEG 

LOCAL  ANGLE  OF  ATTACK  T 027ODEG&PT . OF  UT-0.4)  - 
RAOtAL  POINT  OF  UT-0.4  ?§270DEG)  -  0.400 


3.1  DEG 
3.1  DEG 


APPENDIX  2  -  USING  PROGRAM  'TRANS' 


Program  'TRANS'  (Ref.  11)  was  created  to  plot  and  tabulate 
parameters  from  a  '.DAT'  file,  with  time  as  the  independent  variable. 
When  using  input  files  created  by  'POLAR'  and  'ENERGY'  (Appendix  4), 
the  independent  variable  is  velocity.  Thus  the  time  limit  command  of 
'TRANS',  TIM, now  refers  to  the  range  and  intervals  of  velocity. 

A2.1  Obtaining  a  Speed-Power  Polar 


The  example  below  shows  how  to  run  'TRANS  in  order  to 
create  a  plot  of  total  power  against  velocity  i.e.  a  speed-power 
polar  (Fig.  3).  The  input  file  used  was  created  using  'POLAR'  in 


RU  PUB: TRANS 


[TRANS  version  date  12-NOU-SOI 

I/P  FILENAME  -  POLAR 

BLADE  ELEMENT  METHOD 
POWER  -  SPEED  POLAR 

I/P  FILE  RECORDED  ON  28-Sep-83  AT  09:09 

INTEGN  INT  -  0 . 0000E+00 ;  RUN  CPU  TIME  -  4.02  SEC. 

TIME  FROM  0 . QOOOE+OO  TO  1.5Q00E+Q2  IN  STEPS  OF  2.0OG0E*QO 
(Now  interpreted  as  velocity  range  from  0  to  150  kts  in  steps  of  2  kts) 
♦SPA 

IS  SPACING  BETWEEN  PLOTS  REQRD  :  N 
♦SCA 

BLK  NO.  -1  DENOTES  INDEP  UARIABLE 

ARE  PLOT  SCALE  LIMITS  TO  BE  REAO  FROM  DSK  :  N 

IS  TTY  LISTING  OF  LIMITS  REQRD  :  N 

ARE  MODIFICATIONS  REQRD  :  Y 

BLK.  LOWER,  'JP°ER 
1,0,1500 

♦PLS 

STRIP  PLOTS  : 

BLKS 
1 

TO  SPECIFY  NO.  OF  X  UN  ITS/ INCH,  TYPE  0  FOR  X 
LENGTH  OF  AXES  IN  INCHES;  X,  Y  -  0,6 
MIN  X,  NO.  OF  X  UNITS/INCH  -  0,30 
ARE  SYMBOLS  REQRD  FOR  PLOTS  :  N 
LINE  KEY  (0  GIUES  DEFAULT)  -  0 
♦GOE 

♦*  RUNNING  ♦* 

♦EX  I 

END  OF  EXECUTION 

CPU  TIME:  3.20  ELAPSED  TIME:  1:7.54 
EXIT 


A2.2  Tabulated  Printout 

The  example  below  shows  how  'TRANS 1  is  run  to  create  a 
file  ' POLAR . COL ' ,  containing  tabulated  data  over  a  range  of  velocities. 


.RU  PUB: TRANS 


[TRANS  version  date  12-NOV-801 
I/P  FILENAME  -  POLAR 


BLAOE  ELEMENT  METHOD 
POWER  -  SPEED  POLAR 


I/P  FILE  RECORDED  ON  28-Sep-83  AT  09:10 

INTEGN  INT  -  0.0000E+00;  RUN  CPU  TIME  -  4.02  SEC. 

TIME  FROM  0 . QQ00E+00  TO  1.50Q0E*02  IN  STEPS  OF  2.0000E+00 
♦TIM 

TIME  PARAMS;  LOWER,  UPPER,  INTERVAL  -  60,80,2 
♦PRC 

PRINTING  IN  COLUMNS  : 


BLKS 

A 


IS  O/P  TO  TTY  REQRD  :  N 
♦GOE 

♦♦  RUNNING  ♦♦ 

♦TIM 

TIME  PARAMS;  LOWER,  UPPER,  INTERVAL  -  100,120,2 
♦PRC 

PRINTING  IN  COLUMNS  : 

BLKS 

A 


IS  O/P  TO  TTY  REQRD  :  N 
•GOE 

♦*  RUNNING  ♦♦ 


♦EX  I 

END  OF  EXECUTION 

CPU  TIME:  3.62  ELAPSED  TIME:  57.10 
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APPENDIX  3 


USING  PROGRAM  ’ATMOS' 


The  atmospheric  density,  required  for  'POLAR* ,  and  the  speed 
of  sound,  required  for  'ENERGY',  may  be  found  using  the  program  'ATMOS'. 
(Reference  12).  For  ARL  users,  'ATMOS'  is  available  on  Dectape  D/039  for 
use  on  the  DEC-10  computer,  and  Magtape  M/1012  for  use  on  the  VAX-11. 

The  following  types  of  atmosphere  may  be  considered  by 
setting  the  KEYAIR  flag  appropriately. 

KEYAIR  =  1  :  ICAO  Standard  Atmosphere 

KEYAIR  =  2  :  ICAO  Sea-level  conditions  at  all  times 

KEYAIR  =■  3  s  Off-Standard  ICAO  atmosphere 

KEYAIR  =  4  s  ARDU  Tropical  atmosphere 

KEYAIR  *  5  :  ARDU  Sea-level  conditions  at  all  times 

KEYAIR  *  6  s  Off-standard  ARDU  atmosphere 

Examples  of  running  'ATMOS'  are  shown  below,  demonstrating 
its  use  for  I5A+15"C  conditions. 


.RU  ATMOS 

SET  ATMOSPHERIC  FLAG,  KEYAIR  (1,2  , 3, 4, 5  OR  6 )  :  3 
SINGLE  CALCULATION,  OR  TABLE  Cl  OR  2):  2 
♦STATE  ALTITUDE  (IN  FEET):  20000 
TEMPERATURE  OF  THE  DAY,  TDAY  (IN  DEG.  C):  30 
QNH  OF  THE  DAY  (IN  MILLIBARS):  1013.25 
HEIGHT  OF  THE  AIRFIELD  REFERENCE  POINT,  HAFR:  0 
PAUSE 

T^pe  C  to  Continue,  X  to  Exit,  T  to  Trace. 

END  OF  EXECUTION 

CPU  TIME:  0.51  ELAPSED  TIME:  1.26 
EXIT 


♦  'ATMOS'  calculates  results  in  intervals  of  1000ft 
up  to  stated  altitude. 


The  results  are  stored  in  file  'ATMOS. OUT 


Column  1  is  altitude  in  feet 
Column  2  is  temperature  in  Kelvin 
Column  3  is  pressure  in  lbs/ft2 
Column  4  is  density  in  slug/ft3 
Column  5  is  speed  of  sound  in  ft/s. 


TV  ATMOS. OUT 
0.0000000 
1000.0000000 
2000.0000000 
3000.0000000 
4000.0000000 
5000.0000000 
6000.0000000 
7000.0000000 
8000.0000000 
9000.0000000 
10000.0000000 
11000.0000000 
12000.0000000 
13000.0000000 
14000.0000000 
15000.0000000 
16000.0000000 
17000. 0000000 
18000.0000000 
19000.0000000 
20000.0000000 


303.1499939 

301.1687927 

299.1875916 

297.2063904 

295.2251892 

293.2439880 

291.2627869 

289.2815857 

287.3003845 

285.3191833 

283.3379822 

281.3567810 

279.3756104 

277.3944092 

275.4132030 

273.4320068 

271.4508057 

269.4696045 

267.4884033 

265.5072021 

263.5261012 


2116 

2040 

1967 

1896 

1827 

1760 

1695 

1632 

1571 

1512 

1455 

1399 

1345 


.2199707 

.8562012 

.6798096 

.6431885 

.6977539 

.7960205 

.8903809 

.9365234 

.8883057 

.7016602 

.3320313 

.7375488 

.3762207 


1293.7042236 

1  O  Q  ~ 

1194 ! 2706299 
1146.9234653 
1101. 1173216 
1056 . 7994385 
1013.9373779 
972.4935654 


0.0022593 
0.0021932 
0.0021285 
0.0020654 
0.0020036 
0.0019433 
0.0018844 
0.0018269 
0.0017707 
0.0017159 
0. 0016624 
0.0016101 
0.0015591 
0 . 0015094 
0 . 0014609 
0  .  0  0 1 4  i  3  6 
0.0013675 
0.0013225 
0. 0012787 
0.0012360 
0 . 0011943 


1145.1402588 
1141.3920898 
1137.6317139 
1133.8587646 
1130.0732422 
1126.2750244 
1122.4638672 
1118.6398926 
1114.8026123 
1110.9522705 
1107. 0883789 
1103.2109375 
1099.3200684 
1095.4151611 
1 04 i . 4a-337° 
1037 . 5=333-5 
1083 .6163207 
1079.6542733 
1075.6732227 
1071.6872559 
1067.6812334 


The  example  below  shows  the  running  of  'ATMOS'  when  a  single 
calculation  is  required: 


a 


I 

i 


.RU  ATMOS 

SET  ATMOSPHERIC  FLAG,  KEVAIR  <1,2, 3, 4, 5  OR  6):  3 

SINGLE  CALCULATION,  OR  TABLE  <1  OR  2):  1 

STATE  ALTITUOE  (IN  FEET):  15500 

TEMPERATURE  OF  THE  DAY,  TDAY  (IN  DEG.  C):  30 

QNH  OF  THE  DAY  (IN  MILLIBARS):  1013.25 

HEIGHT  OF  THE  AIRFIELD  REFERENCE  POINT,  HAFR:  0 

PAUSE 

T^pe  C  to  Continue,  X  to  Exit,  T  to  Trace. 


END  OF  EXECUTION 

CPU  TIME:  0.37  ELAPSED  TIME:  1.04 
EXIT 


TY  ATMOS. OUT 
15500.0000000 


272.4414063  1170.4053955 


0.0013904 


1085.5915527 


APPENDIX  4  -  PROGRAM  'ENERGY* 


As  mentioned  in  the  introduction,  program  'ENERGY*  uses  the 
same  relatively  simple  energy  balance  equations  as  'MCEP'  to  estimate  the 
power  components  in  steady  level  flight.  However,  the  profile  power 
equation  used  in  'MCEP'  applies  only  to  two-bladed  rotors  with  zero  twist 
This  equation  has  been  modified  to  apply  to  helicopters  having  any  number 
of  rotor  blades,  but  still  with  zero  twist. 

MCEP 


ENERGY 


where  C  _  Po  _  T _ 

o  pffR2  (QR)  3  an  CT  “  pirR2  (£2R)  2 

It  can  be  seen  that  the  above  equations  are  identical  for 
B  =  0.95,  which  is  a  typical  value  for  a  two-bladed  rotor  (Ref.  13). 
'ENERGY'  calculates  tip  loss  factor  'B'  using  the  standard  equation: 


In  addition  to  the  'MCEP*  relationships,  equation  24  was 
incorporated  into  'ENERGY'  to  enable  estimation  of  instantaneous  rate 
of  climb  for  a  given  shaft  power.  Capability  to  estimate  range  and 
endurance  has  also  been  included  using  the  Breguet  equations. 

(Equations  25  and  26) . 

When  a  speed-power  polar  is  required  'ENERGY'  outputs  a  file 
'ENRG.DAT'  which  is  compatible  with  'TRANS'  (Ref.  11).  If  range, 
endurance  or  one  specified  velocity  is  required  'ENERGY'  outputs  a  file 
'ENERGY. OUT' .  'ENERGY'  also  has  a  file  renaming  feature  similar  to 
'POLAR'  (Appendix  1).  Here  however,  files  are  renamed  ' ENERG1 . OUT '  or 
'ENRG1.DAT'.  Up  to  five  consecutive  runs  may  be  made  without  overwriting 
previous  files,  (as  with  'POLAR') 


-r-  <1  ♦«.«.■)&  ♦*!  fe  £»♦«,£.  £n 


-r  -  «  ♦  4-6“J>t|2  +  5i  43?  <r>  *  w 


A4.1  Loading  'ENERGY'  on  the  DEC-10 


For  ARL  users,  a  description  is  shown  below  of  how  to  load 
program  'ENERGY'  on  the  DEC-10  computer.  'ENERGY'  is  also  available  on 
the  VAX-11  computer  and  is  loaded  in  a  similar  fashion. 


.PLE  MX  D/039 

MESSAGE  TRANSMITTED  at  14:19 


"j  JOPR:  -  D/039  MX  DTA1 
.AS  DTA1 
DTAOOl  assigned 


.DIR  DTA1 : 


Tape  ID:  D/039 
Free:  428  blks, 
POLAR  . CCL  1 

POLAR  .FOR  60 

DING  .REL  1 

ENERGY . CCL  1 
ENERGY. FOR  32 
SLADE  . FOR  43 
T5UB81 • REL  10 

rsi  i  3-1 

-5tygs  .FCR  15 


13  files 
20-Ap  r-83 
7-Sep-33 
14-Apr-83 
20-Apr-33 
7-Sep-33 
6- Jun-33 
30-Mar-33 
a.  sjov-  7-5 
25- Jun-33 


•COPY/X-OTA1 : ENERGY. *,TSUS31 . * .CPU. * , DING. * 

.COM  ENERGY 
FORTRAN:  ENERGY 
ENERGY 
CONTRL 
NEW 
POUJ 
CLIMB 
ENDRNG 

. TY  ENERGY. CCL 
ENERGY , TSU88 1 , CPU .DING 
ENERGY/SAU- 
/GO 

.R  LINK 
*@ENERGY 
EXIT 

.DIR  ENERGY . EXE 

ENERGY  EXE  44  <057>  7-Sep-83 


DSKD:  C1031, 1063  ] 


When  a  speed-power  polar  is  required,  a  file  'ENERGY. OAT' 
is  created  by  'ENERGY'  which  is  suitable  for  use  with  'TRANS'. 

An  example  of  running  'ENERGY'  is  demonstrated  below, 
together  with  the  subsequent  running  of  'TRANS'. 

For  this  case,  the  drag  divergent  Mach  Number  is  assumed  to 
be  0.75,  and  the  thrust  divergence  coefficients  are  0.1  and  0.2.  (TCI 
and  TC2  respectively) .  Definitions  of  induced  velocity  factor,  drag 
divergent  Mach  Number,  TCI  and  TC2  may  be  found  in  'MCEP'  documentation. 
(Refs.  3  and  4]. 

.RU  ENERGY 

TITLE  (TWO  LINES  OF  UP  TO  60  CHARACTERS) 

:  ENERGY  METHOD  (MCEP) 

SPEED  POWER  POLAR 

ARE  UNITS  IN  IMPERIAL  OR  METRIC  ( I  OR  M)  ?  M 

ALL  UP  WEIGHT  (N  OR  LS)  -  ?  42718 

ROTOR  TIP  SPEED  (M/S  OR  FT/S)  -  ?  219.6 

ROTOR  RADIUS  (M  OR  FT)  -  ?  6 . 4 

ROTOR  SLACE  CHORD  CM  OR  FT)  -  ?  3 . 395 

NUMBER  OF  ROTOR  BLAOES  -  ?  4 

2D  LIFT  CURUE  SLOPE  (/RAD)  -  ?  5.73 

INOUCED  VELOCITY  FACTOR  -  ?  2.14 

SPEED  OF  SOUNO  (M/S  OR  FT/S)  -  ?  349 

EQUIVALENT  FLAT  PLATE  AREA  (M**2  OR  FT**2 )  -  ?  1.58 

DRAG  DIVERGENT  MACH  NUMBER  -  ?  0.75 

DRAG  POLAR  CO-EFFICIENT  (DELO)  -  ?  0.0087 

DELTA1  (/RAO)  -  ?  -0.0216 

DELTA2  (/RAD**2)  -  ?  0.4 

AIR  OENSITY  (KG/M**3  OR  SLUG/FT**3 )  -  ?  1.1644 
THRUST  DIVERGENCE  CO- EFFICIENTS  (TCI)  -  ?  0.1 
TC2  -  ?  0.2 

IS  RANGE/ENDURANCE  REQUIRED  ?  (Y  OR  N)  :  N 
IS  A  SPEED-POWER  POLAR  REQUIRED  ?  (Y  OR  N)  :  Y 
MAXIMUM  VELOCITY  (KNOTS)  -  ?  150 
VELOCITY  INTERVAL  (KNOTS)  -  ?  2 

END  OF  EXECUTION 

CPU  TIME:  1.78  ELAPSED  TIME:  1:6.76 
EXIT 


'TRANS'  has  a  plotting  facility  whereby  a  previously  created 
file  ' POLAR . DAT '  (using  'POLAR')  can  be  used  with  the  resulting  'ENRG.DAT' 
file  to  produce  an  overlay  plot.  This  plot  (Fig.  4)  compares  speed-power 
polars  obtained  using  the  two  methods. 

.RU  PUB: TRANS 

CTRANS  version  date  12-NOU-80] 

UP  FILENAME  -  POLAR 

BLADE  ELEMENT  METHOD 
POWER  -  SPEED  POLAR 

UP  FILE  RECORDED  ON  27-0ct-83  AT  17:11 

INTEGN  INT  -  O.OOOOE+OO;  RUN  CPU  TIME  -  4.02  SEC. 

TIME  FROM  0 . 0000E+0Q  TO  1.5000E+02  IN  STEPS  OF  2.0000E*00 
♦SPA 

IS  SPACING  BETWEEN  PLOTS  REQRD  :  N 
♦SCA 

8LK  NO.  -1  DENOTES  INDEP  UAR I ABLE 

ARE  PLOT  SCALE  LIMITS  TO  BE  REhO  FROM  DSK  :  N 

IS  TTY  LISTING  OF  LIMITS  REQRD  :  N 

ARE  MODIFICATIONS  REORO  :  Y 

5LK  ,  uCWER  ,  .JPPES 
1,0,1500 

♦PLS 

STRIP  PLOTS  : 

8LKS 

1 

TO  SPECIFY  NO.  OF  X  UN  ITS/ INCH,  TYPE  0  FOR  X 

LENGTH  OF  AXES  IN  INCHES;  X,  Y  -  0,6 

MIN  X,  NO.  OF  X  UNITS/INCH  -  0,30 

ARE  SYMBOLS  REQRD  FOR  PLOTS  :  N 

LINE  KEY  (0  GIUES  DEFAULT)  -  0 

♦GOE 

**  RUNNING  *♦ 

♦REP 

UP  FILENAME  -  ENRGY 
ARE  SYMBOLS  REQRD  FOR  PLOTS  :  N 
LINE  KEY  (0  GIUES  DEFAULT)  -  0 
♦♦  RUNNING  ** 

UP  FILENAME  - 
END  OF  EXECUTION 

CPU  TIME:  5.72  ELAPSED  TIME:  2:26.10 


3 


Running  'ENERGY'  for  one  particular  airspeed 

If  desired,  'ENERGY'  may  be  run  for  one  particular  airspeed, 
as  shown  below.  The  results  apply  to  level  flight,  the  climb  rate  being 
the  instantaneous  climb  rate  possible  appropriate  to  excess  power,  not 
the  sustained  climb  rate.  The  results  for  this  operation  are  output  on 
to  file  ' ENERGY . OUT ' . 


.RU  ENERGY 

TITLE  (TWO  LINES  OF  UP  TO  60  CHARACTERS) 
i  ENERGY  METHOD  (MCEP) 

RATE  OF  CLIMB 

ARE  UNITS  IN  IMPERIAL  OR  METRIC  (I  OR  M)  ?  M 

ALL  UP  WEIGHT  (N  OR  LB)  -  ?  42718 

ROTOR  TIP  SPEED  (M/S  OR  FT/S)  -  ?  218.6 

ROTOR  RADIUS  CM  OR  FT)  -  ?  6.4 

ROTOR  BLAOE  CHORD  (M  OR  FT)  -  ?  0.395 

NUMBER  OF  ROTOR  BLADES  -  ?  4 

2D  LIFT  CURUE  SLOPE  (/RAO)  -  ?  5.73 

I.NCUCZD  /SLCCtTY  FACTOR  ■  '?  2.14 

SPEED  OF  SOUND  (M/S  OR  FT/S)  »  ?  349 

EQUIUALENT  FLAT  PLATE  AREA  (H**2  OR  FT**2)  -  ?  1.53 

DRAG  DIUERGENT  MACH  NUM8ER  -  ?  0.75 

DRAG  POLAR  CO-EFFICIENT  (DELO)  -  ?  0.0087 

DELTA 1  (/RAD)  -  ?  -0.0216 

DELTA2  (/RAD**2)  -  ?  0.4 

AIR  DENSITY  (KG/M**3  OR  SLUG/FT*+3 )  -  ?  1.1644 
THRUST  DIUERGENCE  CO-EFFICIENTS  (TCI)  -  ?  0.1 
TC2  -  ?  0.2 

IS  RANGE/ENDURANCE  REQUIRED  ?  (Y  OR  N)  :  N 
IS  A  SPEED-POWER  POLAR  REQUIRED  ?  (Y  OR  N)  :  N 
UELOCITY  IN  KNOTS  (0  KNOTS  IS  MINIMUM  ALLOWABLE)  ?  :  50 
WHAT  IS  POWER  AUAILABLE  AT  MAIN  ROTOR  SHAFT  (HP)  ?  :  1275 

STOP 

ENO  OF  EXECUTION 

CPU  TIME:  1.17  ELAPSED  TIME:  42.26 
EXIT 


ENERGY  METHOD  CMCEP) 

RATE  OF  CLIMB 

ALL  UP  WEIGHT  -  42713.0000000  N 

AIR  DENSITY  -  1.1644000  KG/M**3 

TIP  SPEED  -  218.6000000  M/'S 

ROTOR  RADIUS  -  6.4000000  M 

ROTOR  CHORD  -  0.3950000  M 

NUMBER  OF  ROTOR  BLADES  -  4.0000000 

LIFT  CURUE  SLOPE  -  5.7300000  /RAD 

BLADE  SOLIDITY  -  0.0785828 

EQUIUALENT  FLAT  PLATE  AREA  -  1.5800000  M**2 

INDUCED  UELOCITY  FACTOR  -  2.1400000 

DRAG  POLAR  CO-EFF I C I ENT ( DELTA0 ) ■  0.0087000 

DELTA  1  -  -0.0216000  /RAD 

DELTA  2-  0.4000000  /RAD**2 

AOUANCE  RATIQ(MU)  -  0.1176670 

FLAT  PLATE  DRAG  -  608.6095000  N 

THRUST  -  42722.3300000  N 

DISC  AREA  -  128.6796000  M**2 

SPEED  OF  SOUND  -  349.0000000  M/S 

DRAG  DIUERGENT  MACH  NO.  -  0.7500000 

TCI  •  0.1000000 

TC2  -  0.2000000 

INOUCED  UELOCITY  -  5.4247720  M/S 

COMPRESS  I S  I L I  TV  POWER  -  0. 0000000. -sP 

INOUCED  POWER  -  325.6911000  HP 

PARASITE  POWER  -  20.9960500  HP 

PROFILE  POWER  -  215.3766000  HP 

TOTAL  POWER  -  562.0638000  HP 

AUAILABLE  SHAFT  POWER  -  1275.0000000  HP 

CLIMB  POWER  ■  712.9362000  HP 

RATE  OF  CLIMB  -  746.6153000  M/M IN 

STALL  POWER  -  0.0000000  HP 

AIRSPEED  -  50.0000000  KT 


Estimating  Range  and  Endurance  using  ’ENERGY' 


To  estimate  the  range  and  endurance  for  a  given  helicopter 
is  run  as  in  the  previous  example ,  with  the  following  changes 


>»wi  nJ"  -C 


IS  RANGE/ENDURANCE  REQUIRED  ?  (Y  OR  N)  :  Y 
UELOCITY  INTERUAL  (KNOTS)  -  ?  2 

SPECIFIC  FUEL  CONSUMPTION  ?  CKG/HP-HR  OR  LB/HP-HR)  :  0.504 
WEIGHT  OF  FUEL  CARRIED  ?  (N  OR  LB)  :  6400 

STOP 

END  OF  EXECUTION 

CPU  TIME:  1.35  ELAPSED  TIME:  48.02 
EX  I T 


The  output  file  is  shown  below: 


,TY  ENERGY. OUT 

Energy  method  cmcep) 

RANGE/ENDURANCE 

ALL  UP  WEIGHT  -  42718.0000000  N 

AIR  DENSITY  -  1.1644000  KG/M**3 

TIP  SPEED  -  218.6000000  M/S 

ROTOR  RADIUS  -  6.4000000  M 

ROTOR  CHORD  -  0.3950000  M 

NUMBER  OF  ROTOR  BLADES  -  4.0000000 

LIFT  CURUE  SLOPE  -  5.7300000  /RAD 

BLAOE  SOLIDITY  -  0.0785828 

EQUIUALENT  FLAT  PLATE  AREA  -  1.5800000  M*+2 

INDUCED  UELOCITY  FACTOR  -  2.1400000 

DRAG  POLAR  CO-crr  ICIENTOELTAO  J-  J.J037000 

OELTA  1  -  -0.0216000  /RAO 

DELTA  2-  0.4000000  /RAD**2 

S.F.C.  -  0.504  KG/HP-HR 

FUEL  WEIGHT  -  6400.0  N 


RESULTS  FOR  ENDURANCE  OF  2.75  HOURS,  ARE  AS  FOLLOWS  : 

AOUANCE  RATIQ(MU)  -  0.1647338 

FLAT  PLATE  DRAG  -  1192.8750000  N 

THRUST  -  42734.6500000  N 

DISC  AREA  -  128.6796000  M**2 

SPEED  QF  SOUNO  -  349.0000000  M/S 

DRAG  DIUERGENT  MACH  NO.  -  0.7500000 

TCI  -  0.1000000 

TC2  -  0.2000000 

INDUCED  UELOCITY  -  3.9283030  M/S 

COMPRESSIBILITY  POWER  -  0.0000000  HP 

INOUCED  POWER  -  225.1538000  HP 

PARASITE  POWER  -  57.6131700  HP 

PROFILE  POWER  -  227.7853000  HP 

TOTAL  POWER  -  510.5522000  HP 

STALL  POWER  -  0.0000000  HP 

AIRSPEED  -  70.0000000  KT 

1  RESULTS  FOR  RANGE  OF  440.8  KM,  ARE  AS  FOLLOWS  : 


AOUANCE  RAT IO (MU)  -  0.2494540 

FLAT  PLATE  ORAG  -  2735.3350000  N 

THRUST  -  42805.4900000  N 

DISC  AREA  -  128.6796000  M**2 

SPEED  OF  SOUND  -  349.0000000  M/S 

DRAG  DIUERGENT  MACH  NO.  -  0.7500000 

TCI  ■  0.1000000 

TC2  -  0.2000000 

INOUCED  UELOCITY  -  2.6054610  M/S 

COMPRESSIBILITY  POWER  -  16.1808900  HP 

INDUCED  POWER  •  149.5816000  HP 

PARASITE  POWER  -  200.0531000  HP 

PROFILE  POWER  •  260.6611000  HP 

TOTAL  POWER  -  626.4767000  HP 

STALL  POWER  -  0.0000000  HP 

AIRSPEED  -  106.0000000  KT 


Figure  1  i  Forces  acting  on  helicopter  in  steady  flight 


Axis  of 
No  Feathering 


Figure  2  i  Hade  Cleaent  Representation 


rigur*  4  *  Comparison  Of  'POLAR'  And  ' CHER6Y'  For  Holicoptor 
Having  Main  Rotor  Blados  With  Zoro  Dtgrtts  Twist 


Fi^ur*  6  •  Comparison  Of  'PQLRR',  Csssow  1  Tapscott  BnO  Tannor 
For  A  Main  Rotor  Having  Zaro  Stfrtu  Twist 
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